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Abstract 


I 


In  this  paper  a  new  modeling  methodology  to  characterize  failure  processes  in  Time-Sharing 
systems  due  to  hardware  transients  and  software  errors  is  summarized.  The  basic  assumption  made  is 
that  the  instantaneous  failure  rate  of  a  system  resource  can  be  approximated  by  a  deterministic 
function  of  time  plus  a  zero-mean  stationary  Gaussian  process,  both  depending  on  the  usage  of  the 


resource  considered.  The  probability  density  function  of  the  time  to  failure  obtained  under  this 
assumption  has  a  decreasing  hazard  function,  partially  explaining  why  other  decreasing  hazard 


function  densities  such  as  the  Weibull  fit  experimental  data  so  well.  Furthermore,  by  considering  the 
Kernel  of  the  Operating  System  as  a  system  resource,  this  methodology  sets  the  basis  for 
independent  methods  of  evaluating  the  contribution  of  software  to  system  unreliability,  and  gives 
some  non  obvious  hints  about  how  system  reliability  could  be  improved.  A  real  system  has  been 
characterized  according  to  this  methodology,  and  an  extremely  good  fit  between  predicted  and 


observed  behavior  has  been  found.  Also,  the  predicted  system  behavior  according  to  this  methology 


is  compared  with  the  predictions  of  other  models  such  as  the  exponential,  Weibull,  and  periodic 


failure  rate. 
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1.  Introduction 


Fault-tolerance  is  the  attribute  which  allows  the  system  to  continue  with  its  program-specific 
behavior  after  the  occurrence  of  a  failure.  As  reported  in  [Avizienis  78],  full  fault-tolerance  is  currently 
found  only  in  a  few  systems  which  have  exceptionally  strict  reliability  or  availability  requirements.  The 
benefits  of  fault-tolerance  are  not  yet  available  to  the  great  majority  of  users.  [Avizienis  78]  proposes 
five  near  term  goals  to  correct  this  situation :  1 

1.  Acceptance  of  fault-tolerance  as  an  attainable  and  desirable  attribute  of  computing 
systems. 

2.  Development  of  cost  and  benefit  measures 

3.  Development  of  fault-tolerance  specifications  and  acceptance  tests 

4.  Creation  of  experimental  systems 

5.  Advanced  research  in  fault-detection,  recovery,  and  modeling  techniques 

The  work  presented  in  this  report  describes  some  new  results  relating  to  points  2,  3,  and  5  above 
and  proposes  further  research  to  be  done  to  completely  characterize  a  new  modeling  methodology. 
This  methodology  gives  quantitative  relationships  between  performance,  workload,  and  (lack  of) 
reliability  for  digital  computing  systems.  Current  methodologies  for  reliability  assessment  may  provide 
good  models  for  explaining  and  predicting  the  behavior  of  systems  in  the  presence  of  hard  (recurrent) 
faults,  but  the  effect  and  charcterization  of  transient  (non  recurrent)  faults  and  software  (either  design 
or  implementation)  errors  is  still  very  elusive.  These  current  reliability  measures,  while  useful  for 
system  design  trade-offs,  do  not  give  individual  users  a  feeling  of  the  impact  of  unreliability  on 
performance  in  general  purpose  system  operating  under  a  variety  of  workloads.  That  is,  there  are  no 
general  methods  for  a  quantitative  assessment  of  the  benefits  of  fault-tolerance.  Also,  there  are  no 
system  independent  methods  of  evaluating  the  relative  importance  of  software  in  the  unreliability  of 
the  total  system.  And  this  problem  is  particularely  important  given  the  constantly  raising  fraction  of 
software  development  in  the  cost  of  computing  systems.  Finally,  if  reliability  is  postulated  to  be 
workload  dependent,  different  uses  of  a  given  system  may  lead  to  different  reliability  predictions. 
Sections  1-2  thru  1-4  investigate  these  issues  in  more  detail,  while  introducing  a  proposed 
methodology.  First,  let  us  examine  some  basic  definitions. 
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1.1  Definitions 


The  following  concepts  need  to  be  precisely  defined : 

Hardware  Fault  Erroneous  state  of  hardware  due  either  to  failures  of  components  or  to  physical 
interference  from  the  environment 

Hardware  Error  Manifestation  of  a  hardware  fault  within  a  program  or  data  structure. 

Permanent  Hardware  Fault 

Hardware  fault  which  is  continuous  and  stable,  reflecting  an  irreversible  physical 
change  in  the  hardware. 

Transient  Hardware  Fault 

Hardware  fault  due  to  temporary  environmental  conditions. 

Software  Fault  Imperfection  in  the  design  or  implementation  of  a  software  module  such  that  upon 

some  timing  or  value  conditions  in  its  input  data  stream  it  fails  to  accomplish  its 
designed  task. 

Software  Error  Manifestation  of  a  software  fault  within  a  program  or  data  structure. 

System  Failure  Manifestation  of  software  or  hardware  errors  that  force  an  entire  computing 

system  to  suspend  its  operation. 

Since  no  repair  takes  place  after  system  failures  due  to  software  faults  or  transient  hardware  faults, 
the  time  of  system  failure  is  essentially  equal  to  the  system  restart  time.  Since  this  report  is  concerned 
solely  in  modelling  hardware  transient  faults  and  software  faults,  the  words  "system  failure"  and 
"system  restart"  will  be  used  interchangeably  to  describe  the  same  event  in  time. 


1.2  Ine  problem  of  Characterizing  System  Reliability 

Fault-tolerance  has  traditionally  been  characterized  by  relatively  simple  functions  based  on  strict 
assumptions.  The  Reliability  function  R(t)  is  defined  as  the  probability  of  uninterrupted  operation  up  to 
time  t  given  that  all  hardware  was  correctly  operating  at  time  t  ■  0.  R(t)  may  be  used  to  characterize 
either  permanent  or  transient  faults.  The  usual  assumption  is  mad^that  the  failure  rate  is  constant 
and,  for  nonredundant  systems  the  reliability  function  becomes  e  ,  where  \  is  is  the  sum  of  the 
failure  rates  of  all  the  components  in  the  system.  A  very  common  quantitative  measure  is  the  Mean 
Time  To  Failure  (MTTF) 


MTTF  « 


L 


00 

R(t)dt 


(1.1) 
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The  popularity  of  the  MTTF  stems  mainly  from  the  fact  that,  for  nonredundant  systems,  it  is  easily 
estimated  by  dividing  the  time  a  system  is  operational  by  the  number  of  failures  reported.  Other 
reliability  indices  used  to  compare  two  systems  A  and  8,  are  the  Reliability  Improvement  factor  (RIF) 
[Anderson  67] 


RIF 


1ra<» 

i-RbW 


and  the  Mission  Time  Improvement  Factor  (MTIF)  [Bouricious  69] 


(12) 


MTIF  =>  when  RA(TA)  *  R8(TB)  *  R^  (1.3) 

tb 

which  are  useful  only  when  the  system  under  study  must  be  available  for  a  predetermined  period  of 
time  T  called  "mission  time". 

The  concept  of  coverage  [Bouricious  69]  is  defined  as  the  conditional  probability  of  successful 
recovery,  given  that  a  fault  has  occurred.  Although  mathematically  attractive,  coverage  has  proven  to 
be  very  difficult  to  estimate  for  real  systems.  Finally,  if  the  Mean  Time  To  Repair  (MTTR)  is  also 
Known,  an  estimate  of  the  system  usefulness  given  by  the  Availability 


A  ,  -  MS— 
MTTF  +  MTTR 


(1.4) 


These  and  other  measures  traditionally  used  to  compare  systems  do  not  take  into  account  the 
performance  of  the  system  whose  reliability  is  being  measured.  Consider  Table  1-1  which  lists  the 
results  obtained  from  seven  different  experiments  whose  goal  was  explicitly  to  gain  experience  on 
systems  reliability.  Data  for  the  first  system  [Yourdon  72],  was  obtained  from  a  summary  of  failure 
statistics  on  a  Borroughs  5500  over  a  15  month  period  starting  in  April  of  1969.  Limited  information 
about  the  cause  of  each  failure  is  available.  For  instance,  one  of  the  categories  includes  system 
failures  due  to  unexpected  I/O  intercepts.  These  failures  are  recorded  whenever  the  software 
responds  to  an  interrupt  signifying  that  some  I/O  action  has  taken  place,  but  discovers  that  it  has  no 
record  of  having  initiated  such  action.  It  is  thus  an  indication  of  some  form  of  hardware  or  software 
error  but  the  particular  cause  for  the  failure  (hardware  or  software)  remains  unknown.  The  data  for 
the  second  system  was  reported  in  [Lynch  75]  and  comes  from  the  first  thirteen  months  of  operation 
of  an  operating  system  called  Chi/OS  for  the  Univac  1 108  developed  by  the  Chi  Corporation  between 
1970  and  1973.  No  explanation  is  given  about  how  such  an  accurate  decomposition  of  failures  due  to 
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hardware  and  software  could  be  obtained.  [Reynolds  75]  reports  data  obtained  from  a  dual  IBM 
370/165  at  Hughes  Aircraft  Company  over  a  period  of  three  years  installed  to  handle  a  mixed  batch 
and  time  sharing  load.  The  forth  system  is  at  the  Stanford  Linear  Accelerator  Center  (SLAC)  where 
the  main  workload  is  processed  aa  multi-stream  background  batch.  The  system  consists  of  a 
foreground  host  (IBM  370/168)  and  two  background  batch  servers  (IBM  370/168  and  IBM  360/91). 
The  architecture  is  designed  to  be  highly  available  and  reconfigurable.  The  CMU-10A  is  an  ECL  POP- 
10  used  in  the  Computer  Science  Department  at  Camegie-Mellor.  University.  The  data  for  the  CRAY-1 
was  reported  in  [Keller  76],  and  the  data  for  the  three  generic  UNIVAC  systems  was  reported  in 
[Siewiorek  80]. 

Table  1-1  gives,  when  available,  a  Mean  Time  to  restart  (MTTS)  value  in  hours  (that  is,  the  Mean 
Time  to  System  Failure),  a  Mean  Number  of  Instructions  to  Restart  (MNIR)  which  is  an  estimate  of  the 
mean  number  of  instructions  executed  from  system  start  up  until  system  failure,  and  the  percentages 
of  system  failures  that  were  caused  by  hardware  faults,  software  faults,  and  whose  cause  could  not  be 
resolved.  The  information  about  execution  rates  needed  to  compute  the  MNIR  value  was  obtained 
from  [Phister79]. 

Obviously,  the  figures  shown  in  Table  1.1.  do  not  carry  much  information.  A  MTTS  figure  alone  does 
not  tell  the  impact  of  unreliability  on  system  use.  Compare  for  example  the  CRAY-1,  [Russell  78],  with 
the  CMUA,  [Bell  78a].  Although  the  CRAY-1  crashes  twice  as  often  as  the  CMUA,  it  can  operate 
continously  at  rates  above  138  Million  Instructions  Per  Second  (MIPS),  while  the  CMUA  operates  at 
1 .2  MIPS.  Hence  the  CMUA  executes  ~1010  instructions  between  crashes  while  the  CRAY-1  executes 
~1012  instructions  between  crashes.  Inconsistencies  like  this  one  suggest  that  reliability  modelling 
and  measuring  should  be  closely  related  with  the  characterization  of  the  performance  of  the  system 
under  study.  Integrated  performance-reliability  models  have  already  started  to  appear  in  the 
literature.  In  [Meyer  79],  a  performance  measure  called  "performability"  gives  the  probability  that  a 
system  performs  at  different  levels  of  "accomplishment”.  In  [Gay  79],  systems  are  modelled  with 
Markov  processes  in  order  to  estimate  the  probability  of  being  in  one  of  several  capacity  states.  This  is 
a  similar  approach  to  the  one  previoulsy  taken  in  [Beaudry  78],  where  the  concept  of  “computation 
reliability”  was  introduced  as  a  measure  which  takes  into  account  the  computation  capacity  of  a 
system  in  each  possible  operational  state.  Finally,  a  Performance/ Availability  model  for  gracefully 
degrading  systems  with  critically  shared  resources  is  given  in  [Chou  80]. 

However,  most  of  the  above  models  have  been  developed  mainly  for  hard  failures,  that  is,  stable 
failures  that  reflect  an  irreversible  physical  change  in  the  hardware.  Unfortunately,  as  it  has  been 
repeatedly  reported  ( [Fuller  78],  [McConnel  79a],  [Morganti  78],  [Siewiorek  78],  [Ohm  79]),  transient 
failures  occur  at  least  an  order  of  magnitude  more  often  than  hard  failures.  A  cost  effective  analysis 
should  then  consider  transients  as  the  main  reason  for  system  unreliability. 
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System 

MTTS  (hours) 

MNIR 

%HW 

%SW 

%  Unknown 

B5600 

14.7 

39.3 

8.1 

52.6 

Chi/06 
(Univac  1108) 

17 

6.7  1010 

46 

55 

- 

dual 

370/165 

8.86 

2.8  1011 

66 

32 

3 

sue 

20.2 

2.3 1011 

73.3 

21.6 

5.1 

CMU-10A 

10 

4.3 1010 

- 

- 

- 

CRAY-1 

4 

1.9 1012 

- 

- 

- 

UNIVAC 

(Large) 

51 

42 

7 

UNIVAC 

(Medium) 

57 

41 

2 

UNIVAC 

(Small) 

88 

9 

3 

Table  1*1:  Reliability  experience  of  several  comercial  systems. 

MTTS  is  the  Mean  Time  to  restart.  MNIR  is  the  Mean  Number  of 
Instructions  to  Restart 

Simultaneoulsy  with  the  developments  described  above,  qualitative  relationships  between  workload 
and  unreliability  have  also  been  noted.  The  results  published  in  [Beaudry  79]  suggest  a  strong 
dependency  between  workload  and  reliability  of  digital  computing  systems.  And  in  the  paper  by 
[Butner  80],  this  dependency  is  stated  explicitly  claiming  that  a  periodic,  workload-dependent  failure 
rate  is  more  appropriate  to  characterise  the  reliability  of  time-sharing  systems  than  the  classical 
constant  failure  rate  model  traditionally  used.  As  reported  in  [Castillo  80],  if  such  a  dependency  is 
taken  into  account  it  is  possible  to  characterize  the  performance  of  digital  computing  systems 
considering  reliability  as  an  inherent  attribute. 


6 


WORKLOAD,  PERFORMANCE,  AND  RELIABILITY  OF  DIGITAL  COMPUTING  SYSTEMS 


1.3  Software  Reliability 

The  problem  of  software  reliability  assessment  is  part  of  the  more  general  area  of  software  quality 
assessment  [Mohanly  73].  Effective  machanisms  for  measuring  software  quality  are  required  due  to 
the  high  cost  of  software  development  and  maintenance.  By  1985  forecasts  indicate  that  over  90%  of 
the  total  computing  dollars  spent  annually  will  be  for  software  [Horowitz  75].  The  development  of 
techniques  for  measuring  software  reliability  has  been  motivated  mainly  by  project  managers  that 
need  both  ways  of  estimating  the  man-power  needed  to  develop  a  software  system  with  a  given  level 
of  performance  and  techniques  to  detect  when  this  level  of  performance  has  been  reached.  However, 
most  software  reliability  models  presented  up  to  date  are  still  far  from  satisfying  these  two  needs  in  a 
general  context. 

Software  reliability  models  can  be  roughly  grouped  in  four  categories.  The  first  category  would 
include  models  formulated  in  the  time  domain.  These  models  attempt  to  relate  software  reliability 
(characterized,  for  instance,  by  a  MTTF  figure  under  typical  workload  conditions)  to  the  number  of 
bugs  present  in  the  software  at  a  given  time  during  its  development.  Typical  of  this  approach  are  the 
models  presented  in  [Shooman  73],  [Musa  75],  and  [Jelinsky  73].  Bug  removal  should  increase  MTTF 
and  correlation  'of  bug  removal  history  with  the  time  evolution  of  the  MTTF  value  may  allow  the 
prediction  of  when  a  given  MTTF  value  will  be  reached.  An  example  of  the  application  of  time  domain 
models  to  the  development  of  a  real-time  system  is  given  in  [Miyamoto  75].  The  main  disadvantages 
of  time  domain  modek  are  that  bug  correction  can  generate  more  bugs,  and  that  software  unreliability 
can  be  due  not  only  to  implementation  errors  (bugs)  but  also  to  design  (specification)  errors. 

Another  approach  to  software  reliability  modeling  is  based  on  studying  the  data  domain.  The  first 
model  of  this  kind  is  described  in  [Nelson  73].  In  principle,  if  sets  of  all  input  data  values  upon  which  a 
computer  program  can  operate  are  identified,  an  estimated  of  the  reliability  of  the  program  can  be 
obtained  by  running  the  program  for  a  subset  of  input  data  values.  A  more  detailed  description  of  data 
domain  techniques  is  given  in  [Thayer  78].  In  the  paper  by  [Schick  78]  the  time  domain  and  data 
domain  models  are  compared.  However,  different  applications  will  tend  to  use  different  subsets  of  all 
possible  input  data  values,  "seeing"  different  reliability  values  for  the  same  software  system.  This  fact 
is  formally  take  into  account  in  [Cheung  80],  where  software  reliability  is  estimated  from  a  Markov 
model  whose  transition  probabilities  depend  on  a  user  profile.  Techniques  for  evaluating  the 
transition  probabilities  for  a  given  profile  are  given  in  [Cheung  75]. 

The  third  category  includes  models  in  which  software  reliability  (and  software  quality  in  general)  is 
postulated  to  obey  certain  laws  [Ferdinand  74],  [Fitzsimmons  78].  Although  such  models  have 
generated  high  amounts  of  interest,  their  genera)  validity  has  never  been  proven  and,  at  most,  they 
only  give  a  figure  for  the  number  of  bugs  present  in  a  program. 
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Finally,  there  have  been  some  attempts  to  characterize  total  system  reliability  (hardware  and 
software)  in  [Costes  78],  modelling  of  fault-tolerant  software  (through  module  duplication)  in  [Hecht 
76],  and  warnings  about  how  not  to  measure  software  reliability  [Littlewood  79]. 

What  ail  the  above  models  have  in  common  is  that  none  of  them  characterizes  system  behavior 
accurately  enough  as  to  give  to  the  user  a  figure  of  guaranteed  level  of  performance  under  general 
workload  conditions.  They  concentrate  in  estimating  number  of  bugs  present  in  a  program  but  do  not 
give  any  accurate  method  to  characterize  and  measure  operational  system  unreliability  due  to 
software.  There  is  a  wide  gap  between  the  varaiables  that  can  be  easily  measured  in  a  running 
system  and  the  number  of  bugs  in  its  operating  system.  However,  a  cost  effective  analysis  should 
precisely  allow  to  evaluate  the  impact  of  software  unreliability  from  variables  easily  accessible  in  an 
operational  system,  without  knowing  the  details  of  how  the  operating  system  has  been  written. 


1.4  Measuring  Reliability  under  typical  and  atypical  conditions 

The  assumption  here  is  that  reliability  is  a  performance  attritbute  in  the  sense  that  a  lack  of 
reliability  increases  the  expected  value  of  the  system  response  time.  If  such  a  relationship  can  be 
derived  in  a  general  context,  policies  and/or  design  parameters  could  be  used  to  optimize  the 
ultimate  system  performance.  In  this  report,  the  approach  taken  is  that  failure  rate  time  variations 
should  closely  follow  workload  time  variations.  Intuitively,  the  dependency  between  workload  and  lack 
of  reliability  can  be  explained  quite  easily.  Assume  that  we  have  a  constant  failure  rate  for  the  primary 
memory  of  a  digital  computing  system  operating  in  a  stable  environment  under  a  time  sharing  policy. 
That  the  transient  failure  rate  in  a  memory  is  constant  is  a  reasonable  assumption.  There  is 
justification  for  thinking  that  certain  complex  devices  might  follow  an  exponential  failure  law  ( [Barlow 
66],  pp  18-22).  The  physical  characteristics  of  the  memory  IC's  do  not  change  with  time  (at  least 
during  the  effective  life  cycle  of  modem  digital  computing  systems).  We  have  to  look  then  for  the 
origin  of  these  transients  either  in  external  sources,  such  as  radiation,  the  presence  of  noise  (possibly 
impulsive)  in  the  power  supply  or  in  the  limitations  of  the  manufacturing  process.  In  fact,  it  has  been 
reported  in  [Geilhofe  79]  that  MOS  memory  devices  exhibit  non  recurring  bit  failures  caused  by  Alpha 
particles  emitted  from  small  amounts  of  radioactive  elements  present  in  1C  packaging  material.  The 
failure  rate  for  this  kind  of  failures  is  of  course  constant.  Assume  now  that  a  transient  memory  failure 
has  higher  probability  of  leading  to  a  system  crash  when  the  central  processor  is  executing  in  Kernel 
mode  than  when  it  is  executing  in  user  mode.  A  memory  failure  when  the  CPU  is  executing  in  user 
mode  may  affect  a  user  process  but  will  not  crash  the  system.  The  system  failure  rate  due  to  transient 
memory  failures  will  then  depend  on  the  ratio  of  the  number  of  memory  references  while  in  Kernel 
mode  to  the  total  number  of  memory  references  per  unit  time.  Since  it  is  a  well  known  fact  that 
operating  system  overhead  increases  with  workload,  the  previous  ratio  will  also  be  a  nondecreasing 
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function  of  the  system  workload,  increasing  in  turn  the  observed  system  failure  rate.  The  result  is  that 
the  observed  system  failure  rate  due  to  transient  memory  failures  should  be  equal  to  the  sum  of  a 
component  following  the  operating  system  ovehead  variations  in  time  (or  indirectly,  workload 
variations  in  time),  plus  a  constant,  workload  independent  component  (even  if  the  system  is  idle,  there 
may  still  be  memory  errors  that  corrupt,  for  instance,  the  clock  interrupt  subroutine). 

Even  if  the  fact  that  a  computing  system  is  not  always  equally  sensitive  to  the  presence  of  hardware 
errors  is  not  considered,  there  are  still  arguments  to  support  the  idea  that  the  apparent  system  failure 
rate  should  depend  on  the  workload.  The  fact  is  that  in  most  computing  systems,  a  component  failure 
will  be  noticed  only  if  the  component  is  "exercised".  A  time  sharing  system  with  no  load,  spending 
most  of  its  time  in  a  wait  state  and  only  a  fraction  of  the  time  executing  the  clock  interrupt  routine  may 
sustain  several  failures  and  still  not  report  any  errors  if  the  minimal  hardware  configuration  required  to 
execute  these  basic  functions  is  not  affected.  It  is  not  maintained  here  that  failures  will  be  caused  by 
increased  utilization  (although  in  some  cases  this  situation  is  certainly  possible)  but  that  they  will  be 
detected  by  an  increase  in  system  utilization.  This  effect  has  also  been  referred  to  as  "error  latency" 
[Shedletsky  73]. 

Analogous  arguments  lead  to  the  expectation  that  the  rate  of  system  failures  due  to  software 
unreliability  will  depend  on  how  much  the  software  is  exercised.  System  software  failures  are  due  to: 
a)  the  (static)  input  data  to  a  progam  module  presents  some  peculiarities  that  the  program  is  not  able 
of  handling  or,  b)  the  software  is  not  capable  of  handling  some  time  dependent  (dynamic)  sequence  in 
the  input  data  stream.  In  the  case  of  a  time  sharing  system,  the  only  software  capable  of  provoking  a 
system  failure  is  the  Kernel  of  the  Operating  System.  This  software  usually  executes  in  a  privileged 
processor  state  and  a  software  error  that  corrupts  some  critical  information  in  the  Kernel  data 
structures  may  lead  to  a  system  failure.  However,  since  nobody  has  any  a  priori  knowledge  of  what 
these  errors  are,  it  is  less  likely  that  the  system  finds  one  of  these  combinations  in  its  input  stream 
under  low  load  (that  is,  small  amounts  of  input  data  to  process  per  unit  time)  than  in  a  high  load 
situation.  Again,  the  observed  system  failure  rate  has  to  depend  on  the  system  load.  Furthermore, 
upon  correct  system  operation,  a  user  program  is  restricted  to  access  any  resource  for  which  it  has 
not  been  given  explicit  permission  by  the  kernel.  Hence,  it  is  not  necessary  to  worry  about  the  effects 
of  user  programs.  Unfortunately,  a  mathematical  characterization  of  these  phenomena  is  not 
available.  Most  of  the  so  called  software  reliability  models  attempt,  at  most,  to  give  a  figure  for  the 
Mean  Time  To  Failure  of  a  software  system  under  some  "typical"  workload  conditions.  As  will  be  seen 
in  the  following  sections,  the  characterization  of  a  "typical"  workload  is  in  itself  an  important  problem. 

One  of  the  more  important  byproducts  of  considering  a  time  varying  failure  rate  in  which  failures 
can  be  due  either  to  hardware  transients  or  software  design  errors  is  that  the  relative  contribution  of 
software  to  system  unreliability  can  be  estimated  directly  from  the  history  of  system  failures.  From  a 
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software  point  of  view,  the  model  presented  here  is  more  in  the  line  of  the  ideas  exposed  in 
[Uttiewood  79]  in  the  sense  that  the  concepts  of  bug  identification  and  elimination  should  be 
separated  from  reliability  measurement  No  one  cares  about  how  many  bugs  remain  in  a  software 
system  If  the  sytem  operates  at  an  acceptable  level  of  performance.  The  modeling  methodology 
presented  in  this  report  does  not  give  any  solution  to  the  problem  of  improving  software  reliability 
(although  it  gives  some  non  trivial  hints  about  how  that  could  be  done)  but  gives  a  method  to 
characterize  the  distribution  of  the  time  to  failure  due  to  software  under  general  workload  conditions. 

The  formal  characterization  of  performance  of  a  digital  computing  system  may  be  very  elusive.  As 
described  in  [Ferrari  75],  there  are  no  known  system-independent  and  workload-independent 
performance  indices  (two  necessary  properties  to  consider  a  measure  a  universal  measure).  But  the 
average  user  is  only  concerned  in  the  elapsed  time  since  a  computation  is  requested  and  the  correct 
result  is  produced.  This  time  will  depend  on  the  load  of  the  system,  the  operating  system  overhead,  on 
the  probability  that  the  system  foils  and  his  particular  task  has  to  be  restarted,  and  of  course,  on  the 
underlying  hardware  configuration.  It  is  then  an  important  problem  to  establish  formal  quantitative 
relationships  between  workload,  performance,  and  reliability. 

In  summary,  this  report  gives  a  solution  to  the  mathematical  characterization  of  the  relationships 
between  workload,  performance,  and  reliability  due  to  transient  failures  and  software  design  errors. 
The  mathematical  analysis  developed  here  can  be  applied  not  only  to  computing  systems,  but  to  any 
complex  systems  in  which  reliability  is  an  important  characteristic  and  for  which  some  knowledge 
about  workload  variations  is  available.  Since  a  large  class  of  these  systems  operate  under  a  quasi- 
periodic  demand  (such  as  public  transportation  systems,  power  distribution  networks,  time  sharing 
and  some  real-time  computing  systems,  etc.),  the  mathematical  characterization  has  been  developed 
first  for  systems  in  which  the  workload  can  be  characterized  by  a  cyclostationary  stochastic  process 
(a  time  varying  stochastic  processes  with  periodic  mean  and  variance). 

In  Section  2  the  formal  assumptions  made  in  the  characterization  of  the  failure  process  of  a  time- 
shared  computer  are  stated  in  detail.  Also,  general  expressions  are  derived  tor  the  Probability 
Distribution  Function,  Reliability  Function,  and  Hazard  Function  of  the  times  to  hardware  failure, 
software  failure,  and  system  failure  when  the  system  overhead  is  described  by  a  cyclostationary 
process  that  can  be  approximated  by  a  periodic  function  plus  a  "corrected"  zero  mean  Gaussian 
process. 


The  results  presented  in  Section  2  are  elaborated  in  Section  3  where  the  failure  process  of  a  real 
system  is  studied  in  detail,  and  exact  expressions  are  given  that  characterize  the  software,  hardware, 
and  system  reliability  for  that  partcular  system.  In  Section  4  these  results  are  compared  with  the 
available  data  regarding  the  reliability  of  the  system  under  consideration  and  with  the  characterization 
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that  would  result  from  more  traditional  models  such  as  constant  failure  rate  (time  to  failure  expondaHy 
distributed),  WeibuK,  and  periodic  failure  rate.  Finally,  in  Section  5,  a  Hat  of  items  to  be  further 
Investigated  is  proposed,  along  with  some  preliminary  conclusions.  Two  mathematical  derivations 
particularly  tedious  have  been  left  to  appendices  so  not  to  distract  the  reader  with  cumbersome 
details  that  are  not  relevant  to  the  ideas  preesntsd  in  tMa  report 


MATHEMATICAL  CHARACTERIZATION 


11 


2.  Mathematical  characterization 


This  Section  gives  the  mathematical  basis  of  a  model  abie  of  predicting  and  calibrating  the 
unreliability  of  digital  computers  due  to  hardware  transients  and  software  errors.  First,  the  necessary 
definitions  are  given  in  Section  2-1 .  The  assumptions  that  the  systems  to  be  modelled  are  assumed  to 
satisfy  are  stated  in  detail  in  Section  2-2.  In  Section  2-3  a  mathematical  skeleton  is  built  based  on 
these  assumptions.  The  result  is  a  general  expression  for  the  Probability  Distribution  Function  (PDF) 
of  the  time  to  system  failure.  Finally,  in  Section  2.3.  the  general  procedure  for  evaluating  the  maximum 
likelihood  estimates  of  the  model  parameters  is  outlined. 

2.1  Definitions 

A  stochastic  process  {x(t,u);  t€T,  «€G}  is  a  family  of  random  variables  all  defined  in  the  same 
probability  space  Q  and  indexed  by  a  real  parameter  t  that  takes  values  in  a  parameter  set  T  called  the 
index  set  of  the  process.  The  indexing  parameter  t  will  represent  time  in  all  the  processes  presented  in 
this  report  and  T  will  always  be  equal  to  the  real  line  R,  that  is,  only  continuous  time  processes  win  be 
considered.  For  each  fixed  t€R,  x(t,«)  as  a  function  of  u  will  be  a  real  valued  random  variable.  For 
each  u€ Q,  x(t,«)  as  a  function  of  t  will  be  a  real  valued  function  of  time  called  a  realization  of  the 
process.  The  set  of  all  these  time  functions  is  called  the  ensemble  of  the  process. 

Definition  1 :  A  counting  process  {N(t,w);  t>tg}  is  a  stochastic  process  having  the  set 
N+  «  (0,1 ,2,. .,oo}  of  nonnegath/e  integers  as  its  state  space. 

For  each  wcG,  N(t,«)  is  a  piecewise-constant  function  of  t  with  jumps  at  tn(w),  t^u),...,^),  the 
values  of  t1  ,...,  tn  depending  on  the  realization  of  the  process.  Counting  processes  are  always 
associated  with  point  processes,  the  value  of  N(t,u)  for  t.^Ktu1  being  the  total  number  of  "points" 
generated  up  to  t,+1.  All  counting  processes  presented  in  this  report  will  be  associated  to  failure 
processes  of  a  given  system,  the  value  of  N(t,w)  for  t,^Kt.+1  being  the  number  of  system  failures 
detected  up  to  t,+r  A  typical  realization  or  sample  function  of  a  counting  process  is  shown  in  fugure 
2-1. 

Definition  2:  A  Poisson  process  is  a  counting  process  {N(t) ;  t<Jtg}  with  the  following 
three  properties : 

1. PrfNty-O]  -  1 

2.  For  tg<3<t ,  the  increment  N(s,t)  *  N(t)-N(s)  is  Poisson  distributed  with  parameter 
A(t)- A(s),  where  A(t)  is  a  nonnegative,  nondecreasing  function  of  L 


3.  {N(t);te>tg}  has  independent  increments. 
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Figure  2*1:  A  oossible  sample  function  of  a  counting  process. 

Property  3  is  the  distinguishing  property.  It  means  that  for  a  Poisson  counting  process,  the  number 
of  points  in  nonovertapping  intervals  are  statistically  independent  random  variables,  no  matter  how 
large  or  small  the  intervals  are  and  no  matter  how  distant  or  close  they  may  be.  The  function  A(t)  in 
property  2  is  termed  the  parameter  function  of  the  process.  If  A(t)  is  an  absolutely  continuous 
function  of  t,  it  can  be  expressed  as 


t 

\(r)  dr 


(2.1) 


where  X(t)  is  a  nonnegative  function  of  t  for  t^tg.  The  function  X(r)  is  termed  the  intensity  function  of 
the  process  N(t).  At  any  time  t£tg,  the  intensity  function  X(r)  is  the  instantaneous  average  rate  at 
which  points  occur.  If  N(t)  is  a  failure  process  X(t)  is  the  failure  rate  of  the  process. 

Definition  3:  A  Poisson  process  is  said  to  be  homogeneous  when  the  intensity  function 
X(t)  is  a  constant  independent  of  time. 

Definition  4:  Whenever  the  intensity,  function  X(t)  is  not  a  constant  but  a  deterministic 
function  of  time,  the  corresponding  Poisson  process  is  said  to  be  inhomogeneous. 

Definition  5:  Let  x(t)  be  a  stochastic  process  that  is  an  "outside"  process  influencing 
the  evolution  of  a  counting  process  {N(t);t£tQ}.  N(t)  is  a  doubly  stochastic  Poisson  process 
with  intensity  process  (X(t,x(t));t^t0)  if  for  almost  every  realization  of  the  process  x(t),  N(t) 
is  a  Poisson  process  with  intensity  process  function  X(t,x(t)). 


The  process  x(t)  carries  the  information  about  how  the  intensity  process  varies,  and  for  this  reason 
will  be  also  called  the  information  process. 
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Definition  6:  A  stationary  process  (in  the  strict  sense)  is  a  stochastic  process 

{x(t),t€T}  with  the  property  that  for  any  positive  integer  k  and  any  points  t, . \  and  h  in 

T,  the  joint  distribution 

(x(tt) . .x(y) 

is  the  same  distribution  of 

{x(t.,  +  h) . .x^  +  k)} 

Intuitively,  a  process  is  stationary  if  it  has  the  same  joint  statistics  regardless  of  where  the  time 
origin  is  set.  Hence,  if  x(t)  is  a  stationary  Gaussian  process,  the  joint  distribution  function  of 

{x(t,  +  h) . .xftfc  +  h)}  is  a  multivariate  Gaussian  distribution  whose  covariance  matrix  is  independent 

of  h. 

Definition  7:  The  Autocorrelation  function  Rxx(t1,t2)  of  a  process  x(t)  is  defined  as 


R«(W  *  EWVXM 

-/  da2 

where  E{..}  stands  for  expected  value  and  pw, .  w.  .(a^aJ  is  the  joint  probability  density 
function  of  x(t,)  and  x(y.  **** 

If  x(t)  is  stationary  and  real,  Rxx(t1,y  depends  only  on  the  time  difference  r  ■  It^y  and 

HJr)  -  E{x(t+T)x(t)) 

Definition  8:  A  stochastic  process  x(t,u)  is  ergodic  in  the  most  general  sense  if  ail  its 
statistics  can  be  determined  from  a  single  realization  xft.Wg)  of  the  process. 

Loosely  speaking,  a  process  is  ergodic  if  time  averages  (the  only  ones  that  can  be  obtained  from  a 
single  realization  of  the  process)  equal  ensemble  averages  (i.e.  expected  values).  Obviously, 
ergodicity  can  be  defined  with  respect  to  certain  parameters  of  the  process.  Only  ergodicity  with 
respect  to  the  autocorrelation  function  will  be  needed  in  this  report,  which  is  defined  as  follows : 

Definition  9:  A  stochastic  function  is  ergodic  with  respect  to  the  autocorrelation 
function  if 

«„W-tToo  £/Tx(t*T)x(t)dt 


If  ergodicity  of  the  autocorrelation  function  is  satisfied,  the  autocorrelation  function  can  be 
estimated  by  computing  the  above  integral  for  a  finite  record  of  a  single  realization  of  the  process  x(t). 

Definition  10:  A  real  valued,  continuous  time  stochastic  process  is  defined  to  be  a 
cyclostationary  process  with  period  T  if  and  only  if 
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1.  E{x(t)}  .  E{x(t  +  T)> 

2.  E{x(t)x(s)}  ■  E{x(t + T)x(s  +  T)}  V  s,t 

that  is,  it  is  a  stochastic  process  with  periodic  mean  and  autocorrelation  functions. 

Definition  1 1:  A  doubly  stochastic  Poisson  process  will  be  said  to  be  a  cyclostationary 
Poisson  process  if  its  information  process  is  cyclostationary. 

In  summary,  and  as  a  short  introduction,  this  report  summarizes  the  results  obtained  by  assuming 
the  failure  processes  of  Time-Sharing  computing  systems  to  be  characterized  by  cyclostationary 
Poisson  processes. 

2.2  Basic  assumptions  made  in  the  characterization  of  failure 
processes 

First,  the  behavior  of  Mures  in  Time-Sharing  computing  systems  will  be  characterized.  Since  the 
occurrence  of  failures  is  random,  a  necessary  requirement  to  understand  the  process  of  how  a  lack  of 
reliability  affects  the  performance  of  a  system  is  to  find  an  expression  for  the  probability  density 
function  of  the  time  to  failure. 

2.2.1  Characterization  of  the  failure  process 

The  approach  taken  has  been  to  assume  that  the  different  subsystems  failure  processes  can  be 
accurately  modeled  by  cyclostationary  Poisson  processes.  Although  it  is  common  in  reliability  theory 
to  assume  that  failure  processes  are  property  modeled  by  Poisson  processes,  one  may  well  wonder 
why  this  assumption  leads  to  good  results.  There  are  at  least  three  reasons  for  characterizing  failure 
processes  with  Poisson  processes. 

First,  the  conditions  for  a  Poisson  process  are  very  likely  to  be  valid  for  many  physical 
environments.  Qualitatively,  these  conditions  can  be  summarized  as  foNowa : 

•Two  failures  cannot  occur  simultaneously. 

•At  any  time,  there  exists  an  instantaneous  failure  rate  at  which  failures  occur  per  unit  time 
and  such  that  the  value  of  this  instantaneous  failure  rate  is  independent  of  the  past 
history  of  the  system. 

•the  number  of  failures  at  start  time  is  zero. 

(see  [Sneyder  75]  for  a  formal  proof  that  the  above  are  sufficient  conditions  for  a  process  to  be 
Poisson).  If  this  "instantaneous”  failure  rate  is  a  constant,  the  above  three  conditions  define  a 
homogeneous  Poisson  process,  for  which  the  interarrival  times  are  independent  and  exponentially 
distributed  random  variables.  If  the  failure  rate  is  a  deterministic  function  of  time,  a  nonhomogenous 
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Figure  2-2:  Average  number  of  blocks  accessed  in  the  file 
system  as  a  function  of  time  of  day. 


Figure  2-3:  Disks  failures  as  a  function  of  time  of  day. 
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Poisson  process  is  defined.  Finally,  if  the  failure  rate  is  another  stochastic  process,  the  above  three 
definitions  define  a  doubly  stochastic  Poisson  process. 

The  second  reason  for  using  a  Poisson  process  is  that  whenever  we  have  a  point  process  that  is  the 
result  of  pooling  the  points  of  many  independent  point  processes  (  whatever  their  characterization 
may  be),  and  the  component  processes  are  sufficiently  sparse,  the  pooled  process  converge  to  a 
Poisson  process  [Cinlar  72].  This  is  certainly  the  case  of  modem  digital  computing  systems.  The 
complexity  of  a  minicomputer  like  the  PDP-i 1  /40  [Bell  78b]  in  a  minimal  configuration  of  64  Kbytes  of 
memory,  dock,  and  a  terminal  interface  is  on  the  order  of  103 1C  packages.  For  an  supercomputer  like 
the  CRAY-1  [Russell  78],  the  complexity  is  on  the  order  of  10s  1C  packages.  The  average  Mean  Time 
To  Failure  (MTTF)  per  component  is  on  the  order  of  106  hours  (~  103  years)  for  hard  failures  [Hodges 
77].  Hence,  the  system  failure  rate  due  to  transients  is  the  superposition  of  ~103  failure  processes, 
the  probability  of  observing  a  failure  of  any  of  the  component  processes  in  a  meaningful  time  interval 
is  very  small  (of  the  order  of  10*4  for  a  month  interval).  The  fact  that  the  superposition  of  sparse  point 
processes  converges  to  a  Poisson  process  guarantees  that,  independently  of  the  characterization  of 
each  of  the  component  processes,  the  system  failure  process  will  be  very  close  to  a  (non  necessarely 
homogeneous)  Poisson  process. 

Finally  and  most  importantly,  even  if  system  characterization  by  means  of  Poisson  processes  is  only 
approximate,  these  processes  are  very  well  understood  and  fairly  complex  mathematical  tools  exist 


Figure  2-4:  Number  of  blocks  accessed  per  unit  time  in  a  file 
system  during  five  consecutive  weekdays  (millions  of  blocks 
accessed  per  5  minutes). 


MATHEMATICAL  CHARACTERIZATION 


18 


WORKLOAD,  PERFORMANCE,  AND  RELIABILITY  OF  DIGITAL  COMPUTING  SYSTEMS 


That  a  doubly  stochastic  Poisson  process  should  be  used  (that  is,  the  failure  rate  being  another 
stochastic  process)  is  a  fact  suggested  by  the  data  presented  in  Figures  2-2  thru  2-6.  Figure  2-4 
shows  the  values  of  the  number  of  blocks  read  and  written  to  the  Tile  system  of  a  time-sharing  digital 
computing  system  during  five  consecutive  weekdays.  There  is  a  clear  (although  nondeterministic) 
periodicity  in  this  data  .Note,  for  instance,  that  there  is  always  a  peak  after  a  new  day  is  started.  This 
peak  is  due  to  the  backup  of  disks  to  magnetic  tape,  which  the  operator  does  daily  after  midnight  If 
we  average  the  data  for  these  five  days  and  plot  the  profile  for  the  average  disks  use  in  one  day  we 
see  that  there  is  a  common  time  varying  pattern  for  ail  days  (Figure  2-2).  If  we  now  examine  the  one 
day  profile  of  disks  failures  detected  during  the  same  period  of  time  in  Figure  2-3  we  note  a 
remarkable  similarity  between  the  two  plots.  Although  different,  the  plots  in  Figures  2-2  and  2-3 
present  the  main  peaks  and  valleys  at  approximately  the  same  time  of  day.  It  seems  that  in  the  long 
run,  after  averaging  over  a  one  day  period  both  the  failure  rate  and  the  system  usage  variables  show 
the  same  temporal  behavior.  If  such  a  dependency  exists  instantaneoulsy,  that  is,  if  the  failure  rate  at 
a  given  time  depends  on  the  system  load  at  that  time,  it  is  clear  that  the  failure  rate  must  be 
characterized  as  a  stochastic  process,  since  the  load  variations  presented  in  Figure  2-4  cannot  be 
considered  deterministic. 

Figures  2-5  and  2-6  show  the  average  fraction  of  time  in  kernel  mode  for  a  Time-Sharing  system  and 
the  number  of  crashes  detected  during  29  days,  both  plots  as  a  function  of  time  of  day.  Again,  there  is 
some  simlilarity  between  the  two  plots.  The  fraction  of  time  in  kernel  mode  for  a  Time  Sharing  system 
during  five  consecutive  days,  shown  in  Figure  2-7  suggests  a  cyclostationary  process. 

Figures  2-2  thru  2-6  should  be  enough  evidence  to  justify  an  experiment  based  on  the  assumption 
that  failure  rate  is  a  stochastic  process.  Let  \(t)  be  the  value  of  the  instantaneous  failure  rate  at  time 
t.  For  a  doubly  stochastic  Poisson  process,  the  probability  density  function  of  the  time  between 
failures  conditioned  to  a  realization  of  the  process  \(t)  is  given  by  [Sneyder  75] 

/X(T)dT 

. u  (2.2) 


2.2.2  Failure  rate  characterization 

Based  on  the  arguments  in  Section  1 ,  it  will  be  assumed  that  the  instantaneous  value  of  the  failure 
rate  for  a  particular  resource  is  a  nondecreasing  function  of  the  "utilization”  of  that  resource.  For 
instance,  more  failures  per  unit  time  will  be  detected  in  a  file  system  when  the  number  of  blocks  read 
and  written  to  the  disks  per  unit  time  is  near  its  maximum  value  than  when  it  is  used  only  occasionally. 
The  fact  that  system  crashes  occur  more  often  in  periods  of  high  load  has  been  noted  in  [Butner  80]. 


Figure  2-7:  Fraction  of  time  in  kernel  mode,  k(t),  during  five 
consecutive  weekdays. 

The  exact  nature  of  the  functions  relating  resource  utilization  and  failure  rates  may  be  complex, 
different  for  each  resource  and  difficult  to  characterize  from  observed  data.  Since  no  previous 
experience  has  been  reported  of  working  under  these  assumptions,  a  cautious  approach  will  be  taken 
and,  as  a  first  step,  only  linear  relationships  will  be  considered.  In  general  then,  the  failure  rate  A(t)  of 
a  particular  resource  whose  use  is  characterized  by  a  function  u(t)  will  be  given  by 


X(t)  «  a  u(t)  +  b  (2.3) 

where  u(t)  will  be  a  function  such  as  the  ones  shown  in  Figures  2*4  or  2-7.  For  instance,  the  failure 
rate  of  a  file  system  Adk(t)  will  be  given  by  . 

Adk(t)  ■  sdk  b(t)  +  c^  (2.4) 


where  b(t)  is  equal  to  the  sum  of  blocks  read  and  written  to  the  file  system  per  unit  time  as  shown  in 
Figure  2-4.  sdk  is  a  sensitivity  coefficient  relating*  disks  usage  to  failure  rate  and  the  offset  term  c^ 
should  take  care  of  any  possible  drift  in  the  relation  between  usage  and  failure  rate. 


20 


WORKLOAD,  PERFORMANCE,  AND  RELIABILITY  OF  DIGITAL  COMPUTING  SYSTEMS 


The  system  failure  rate,  that  is.  the  rate  at  which  the  system  crashes  and  has  to  be  restarted  from 
scratch  is  not  so  obviously  characterized.  The  protection  mechanisms  provided  by  the  state  of  the  art 
operating  systems  and  computer  architectures  try  to  maintain  continued  system  operation  regardless 
of  individual  component  or  subsystem  failures.  The  fact  is  that  in  most  computers  the  CPU  executes  in 
one  of  several  processing  modes,  each  of  the  modes  having  different  privileges  respect  to  the  overall 
system  control.  A  system  crash  due  to  a  hardware  transient  is  only  possible  when  it  affects  the 
operation  of  the  system  in  the  most  privileged  mode,  the  only  one  able  of  halting  the  entire  system  or 
entering  into  an  infinite  loop  with  no  other  entity  capable  of  correcting  the  situation.  This  most 
privileged  mode  of  operation  is  usually  referred  to  as  the  kernel,  and  the  system  failure  rate  should  be 
a  nondecreasing  function  of  the  fraction  of  time  that  the  system  operates  in  kernel  mode,  that  is, 


+  ®hw 


(2.5) 


where  X^ft)  is  the  system  failure  rate  due  to  hardware  transients,  shw  is  a  sensitivity  coefficient,  k(t)  is 
the  instantaneous  value  of  the  fraction  of  time  that  the  system  operates  in  kernel  mode  and  chw  is  a 
residual,  workload  independent,  failure  rate  (even  if  the  kemei  is  only  slightly  exercised  there  is  the 
possibility  that  a  transient  in  the  main  mamory  will  corrupt  parts  of  the  kernel  data  structures). 

The  system  failure  rate  due  to  software  errors  will  also  depend  on  the  fraction  of  time  that  the 
system  operates  in  kernel  mode  because  the  kernel  of  the  operating  system  is  the  only  software 
capable  of  leading  to  a  system  crash.  However,  when  the  workload  is  very  low,  and  the  kemei 
executes  only  relatively  simple  operations  it  is  to  be  expected  that  this  part  of  the  kernel  will  be  well 
debugged  such  that  the  system  failure  rate  is  zero  for  low  values  of  k(t).  It  will  then  be  assumed  that 
the  software  failure  rate  will  be  zero  for  values  of  k(t)  below  a  threshold  value  kg  and  increase  with  k(t) 
above  kg.  Again,  the  relationship  between  k(t)  and  failure  rate  will  be  assumed  linear  such  that 


X.W  -  O 


othtnxMs 


(2.6) 


where  \Sb(t)  is  the  system  failure  rate  due  to  software  errors,  s^  is  its  sensitivity  coefficient,  and  kg  is 
the  value  of  k(i)  below  wich  this  failure  rate  is  zero.  XM(t)  can  be  rewritten  as 
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(2.7) 


where 


(2.8) 


The  following  expression  can  then  be  obtained  for  the  system  failure  rate,  that  is,  the  rate  at  which 
the  system  crashes  due  either  to  hardware  transients  or  software  errors 


Vf)  3  [shw  +  ssJkW  +  Chw-s^W-SawNi 


(2.9) 


Only  these  three  cases  (the  failure  process  of  a  file  system,  the  system  failure  process  due  to 
transients,  and  the  system  failure  process  due  to  software  errors)  will  be  studied  in  this  report 
Expressions  for  the  probability  density  function  of  the  time  to  failure  and  reliability  function  for  the 
three  cases  are  given  in  Sections  2-2  thur  2-4. 


2.2.3  Workload  characterization 

Something  more  can  be  said  about  the  "utilization"  functions.  Although  being  nonstationary 
processes,  it  is  obvious  that  due  to  the  operational  policies  that  regulate  the  use  of  Time-Sharing 
systems,  they  will  have  a  periodic  behavior.  The  second  hypothesis  that  we  make  is  that  workload, 
and  hence  system  usage  for  time  sharing  systems  can  be  modeled  as  a  cyclostationary  process 
[Gardner  75],  [Gardner  78].  A  cyclostationary  process  is  defined  as  a  second  order  process  with 
periodic  mean  and  autocorrelation  function.  The  periodicity  of  the  mean  is  obvious  from  Figure  2-4 
and  in  fact  it  is  possible  to  make  the  simplifying  assumption  that  the  workload  causing  such  overhead 
can  be  described  by  a  periodic  (hence  deterministic)  function  of  time.  This  is  the  approach  taken  in 
[Butner  80],  where  it  is  expected  that  a  periodic  failure  rate  Poisson  process  will  lead  to  a  more 
accurate  failure  process  characterization  than  a  homomgeneous  Poisson  process  model  (time  to 
failure  exponentially  distributed).  Here,  the  instantaneous  value  of  the  failure  rate  will  be  considered  a 
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random  variable  with  periodic  mean,  and  the  failure  rate  will  be  a  cyclostationary  process.  The  third 
hypothesis  is  that  u(t)  (the  usage  function  of  a  particular  system  resource)  can  be  properly  modeled 
by  adding  a  deterministic,  period  function  of  time  m(t)  plus  a  stationary,  zero  mean,  Gaussian  process. 
That  is, 


u(t)  -  m(t)  +  z(t) 


(2.10) 


such  that  in  general 


X(t)  a  am(t)  +  az(t)  +  b 


(2.11) 


where  m(t)  is  a  periodic,  deterministic  function  of  time  and  z(t)  is  a  stationary,  zero  mean,  Gaussian 
process,  independent  of  m(t).  This  third  hypothesis,  although  attractive,  cannot  be  correct  If  z(t)  is  a 
purely  Gaussian  process,  there  is  a  non-zero  probability  that  X(t)<0  and  the  above  expression  cannot 
be  used  as  a  failure  rate  of  a  Poisson  process.  To  avoid  this  problem  let 


2<t)  «  { 


+  '  umin 


otherwise 


and set 


(2.12) 


u(t)  -  m(t)  +  z(t)-z(t)  (2.13) 

from  where  we  obtain 

X,(t)  -  a{  m(t)  +  z(t)  -  f(t)  ]  +  b  (2.14) 

■  a  m(t)  +  b  +  a(z(t)  -  2(t)]  (2.15) 


In  summary,  the  three  hypothesis  in  which  this  work  is  based  are : 

1.  The  failure  process  of  a  digital  computing  system  can  be  described  by  a  doubly 
stochastic  Poisson  process. 
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2.  The  failure  rate  is  a  linear  function  of  the  operating  system  overhead  (so,  indirectly, 
depends  on  the  system  workload) 

3.  For  computing  systems  with  cyclostationary  workload,  system  overhead  time  variations 
can  be  modeled  as  a  periodic,  deterministic  function  of  time  plus  a  stationary,  zero  mean 
Gaussian  process,  independent  of  the  underlying  periodic  function  and  adequately 
corrected  in  order  to  have  a  positive  failure  rate. 


Note  that  assumption  one  is  much  less  restrictive  than  the  usual  assumption  of  considering  the  time 
between  failures  being  exponentially  distributed  (i.e.,  the  failure  process  is  usually  considered  a 
homogeneous  Poisson  process).  In  later  Sections,  the  insight  gained  in  understanding  system 
behavior  from  dropping  this  oversimplification  will  be  discussed.  Also,  the  implications  of  considering 
(or  not  considering)  assumptions  2  and  3  will  be  discussed. 


2.3  Characterization  of  a  file  system  failure  process 


As  a  first  application  of  the  hypothesis  described  above,  the  failure  process  of  a  file  system  under 
cyclostationary  workload  will  be  studied  in  detail.  The  hypothesis  are  that  the  subsystem  failure  rate  is 
given  by 


\*«  "  8dk  ***>  +  c* 


(2.ie) 


where 


b(t)  -  m^t)  ♦  zdk(t)  -  f^t) 


(2.17) 


The  pdf  of  the  time  to  failure  conditioned  to  a  realization  of  the  process  Xdk(t)  is  given  by 


Pdk(t|Adk(T),to£T£t)  -  A^We' 


\*<T»T 


(218) 


The  general  pdf  is  given  by 


p.«>  •  E0*«> 


,/  v* 


} 


«-19) 
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where  ttw  expectation  is  taken  over  the  ensemble  realizations  of  the  process  A^t).  It  is  shown  in  the 
Appendix  I  that,  under  the  assumption  that  A^t)  is  given  by  (2.14),  the  above  expectation  is  equal  to, 

o«M  •-£{ **<« ' '**'’’* * ' '’***■““'* }  «*» 

The  meanings  and  values  of  each  of  the  parameters  on  which  p(t)  depends  are  described  in  detail  in 
Appendix  I,  and  only  a  summary  wi#  be  given  here,  ^(t)  is  a  periodic  function  of  time,  depending  on 
the  periodic  component  of  A^ft).  The  first  term  in  the  exponent  is  the  variance  of  the  integral  of  z^t), 
and  depends  on  the  autocorrelation  function  of  z^t),  Rn(t).  The  last  term  depends  on  the  mean 
value  of  the  deterministic  part  of  X^ft)  and  the  correction  factor  pdk(s4k,bm|n)  takes  care  of  the 
contribution  of  f^t).  Finally,  it  should  be  noted  that  this  expression  is  only  valid  when  the  second 
derivative  of  the  autocorrelation  function  of  z^t)  at  the  origin  is  finite. 

The  following  expression  can  be  obtained  for  the  Probability  Distribution  Function  of  the  time 
between  errors : 


Pdi»(t<T>  *  f  P»<* 

2  JZJQL  t2  JzfilL  .  [a dkmdh  cdk  •  Prti,(*fiv^,lwn)lT 

-  *dk(0)  e  dk  2  -*dk(T)e  dk  2  dK  dk  dk  dk  * 

“  1  -  ^dk(  t)  2*'  *  I#dkmdk  *  f>dk<sdk,bmm),T 


(2.21) 

(2.22) 

(2.23) 


To  compare  our  model  with  a  real  system  we  still  need  to  estimate  the  parameters  s^  and  c^  from 
observed  data  and  obtain  analytical  expressions  for  the  autocorrelation  function  Ra(r)  and  the 
variance  o2(t)  in  equation  (2.23).  The  general  problem  of  parameter  estimation  for  doubly  stochastic 
Poisson  processes  is  described  in  Section  2*5,  and  a  numerical  procedure  for  estimating  sdk  is  given 
in  section . 


2.4  The  system  failure  process 

The  expression  for  the  system  failure  rate  due  to  hardware  transients  and  software  errors  has  been 
given  in  (2.9),  where  k(t)  is  the  fraction  of  time  that  the  system  operates  in  kernel  mode.  With  the 
hypothesis  that 
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k(t)  -  nyt)  +  Zsy(t)  •  yt)  (2.24) 

(2.9)  can  be  rewritten  as 

\yW  "  [s»w  +  ShJl%W  +  zw(t) '  V«  +  °hw  *  8wrfwW  *  s-ko  <225) 


where 


♦  KO  +  k,,  >40<V%*i 

otharwia* 


An  additional  assumption  has  been  made  here,  that  k(njn<kQ<m|njn.  That  is,  it  is  assumed  that  the  value 
at  which  the  failure  rate  due  to  software  failures  starts  being  nonzero  (kg)  lies  somewhere  between  the 
minimum  value  of  the  periodic  component  of  k(t)  (m^)  and  the  minimum  value  of  k(t)  (k^.  The 
reason  for  this  assumption  is  that  only  in  this  case  a  closed  form  expression  can  be  found  for  the  pdf 
of  the  time  to  system  failure.  Whether  this  assumption  holds  or  not  in  a  real  system  is  checked  later  in 
the  report. 


Again,  the  pdf  of  the  time  to  system  failure  is  given  by 


P«y<t> 


VT)dT 


(2.26) 


Using  the  results  of  Appendix  I,  the  following  expression  is  obtained 

a  r  .  ...  J 

p. 


g 

,  (t)..-2-{*  (t,  e(**w  *  K*sw  +  *hw)may  +  chw‘Psy(8sw  +  shw'kmin)'Psy(9iw'Mt  }  (227) 

•y  9t 


and  the  following  expression  is  obtained  for  the  PDF  of  the  time  to  system  failure 
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P^Kt)  -  1-^(r)e 


2 

<8»w + ®hw^  j  ^  *  *hw^n  sy  ®hw '  ^sy^ssw + “hw^min^  *  gfl) 


Again,  to  completely  characterize  the  failure  process  of  a  real  system,  the  values  of  3W,  shw,  hw,  kg 
need  to  be  estimated  from  the  history  of  failures  of  the  system. 


2.5  Parameter  estimation 

The  general  problem  of  parameter  estimation  for  doubly  stochastic  Poisson  processes  can  be 
stated  as  follows.  Let  {N(t)^>tQ}  be  a  doubly  stochastic  Poisson  counting  process  with  intensity 

\(t,z(t),x),  where  z(t)  is  an  stochastic  process  and  *  •  (xvx2, . ,xm)  is  a  vector  of  unknown 

parameters.  The  occurrence  density  function  that  a  given  realization  of  the  process  has  a  failure  at 
time  tf  if  it  has  been  started  at  time  ]&  is,  given  by 

p(tfptz(T),«T<tf)  .  \(tf,z(tf),i)e\{  X(Tll(T)^dT  (2.29) 

If  we  observe  n  failures  at  times  tf1 . tfn  with  associated  starting  times  ts1 . tsn,  the  probability 

density  function  of  observing  such  set  of  events  is 

P<n,(tf, . tfnir,z(T),ts,<T<tf1,Vi)  -  nr. , P(ts,) X<tf,,z<tf,),*) e^  X(T,Z(T)^dT  (2.30) 

where  P(tSj)  is  the  a  priory  probability  that  the  system  is  started  at  time  ts, .  Taking  the  expectation  with 
respect  the  statistics  of  z(t)  we  can  obtain, 

. tfn|x,ts, . ts„)  -  E  { II", ,  P(ts,)  X(tf|,z(tf,),x) e^  X(T*T)**dT  }  (2.31) 

The  maximum  likelihood  estimate  t'  •  (xj,x2 . .xj^)  of  of  x  in  terms  of  a  particular  realization  of 

the  process  is  by  definition  the  value  of  x  that  maximizes  the  above  density  function  [Meisa  78].  That 

is,  p(n)(tf1 . tfnpf,t8,,i » 1, . ,n)  will  be  maximum  for  it  ■  x\  In  the  cases  presented  in  this  report, 

closed  form  expressions  have  been  obtained  for  the  pdf  of  the  time  to  failure.  They  are  all  of  the  form 
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p(t)  -  H(tf.x)  e 


•H(tf,e,if 


(2.32) 


the  function  to  be  maximized  is  then 

p(n)(tf1 . ty  .  nr.  1  P(tSj)  h(ttp.x)  (2.33) 

Note  that  this  problem  is  equivalent  to  minimazing  the  function 

l(x)  -  X]",  1  H(tfj,tSj,5f)  -  E"s ,  ln(h(tf,.x )]  (2.34) 

subject  to  the  constraints 


h(tf,,?)>0  i»1 . n 


(2.35) 


Since  doeed  form  expressions  for  the  components  of  x  at  the  minimum  are  not  generaly  available, 
this  is  a  typical  nonlinear  programming  problem,  subject  to  nonlinear  inequality  constraints.  Since  this 
problem  will  have  to  be  solved  every  time  that  the  failure  process  of  a  resource  has  to  be  modeled  for 
a  real  system,  particular  care  has  been  taken  in  finding  an  efficient  procedure  for  the  location  of 
minimums  of  functions  of  the  type  (2.34).  In  Appendix  II,  this  procedure  is  described,  along  with 
detailed  descriptions  of  all  the  functions  for  which  it  has  been  used  in  the  evaluation  of  maximum 
likelihood  parameters. 

2.6  The  implications  of  a  workload  dependent  model  in 
software  reliability  evaluation 

A  general  methodology  for  characterizing  system  reliability  in  terms  of  resource  utilization  functions 
has  been  presented  in  the  previous  sections.  First,  the  "typical"  measuring  conditions  have  been 
generalized  to  a  situation  in  which  workload  patterns  are  mapped  into  resource  utilization  functions, 
modeled  by  cyclostationary  processes.  Second,  by  considering  the  kernel  of  the  operating  system  as 
a  system  resource,  an  integrated  hardware/software  reliability  model  has  been  built.  The  assumption 
is  that  the  system  failure  rates  due  to  hardware  transients  and  software  errors  depend  on  the  kernel 
utilization  process.  Third,  since  the  functional  dependencies  of  the  failure  rates  due  to  hardware 
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transients  and  software  errors  with  respect  to  kernel  utilization  are  different,  it  is  in  principle  possible 
to  evaluate  the  relative  contribution  of  each  failure  rate  to  the  unreliability  of  the  total  system. 

Once  the  general  functional  dependency  between  failure  rate  and  kernel  utilization  has  been 
established,  all  it  is  needed  to  completely  characterize  a  real  system  is  to  evaluate  the  maximum 
likelihood  values  of  the  function  parameters.  But  all  it  is  needed  to  evaluate  the  maximum  likelihood 
values  of  these  parameters  is  a  history  of  system  failures.  Hence,  the  contribution  of  software  to 
system  unreliability  can  be  evaluated  just  knowing  the  times  of  a  set  of  system  failures,  without 
needing  any  information  about  how  the  kernel  has  been  written,  let  alone  how  many  bugs  remain  in  it 
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3.  Failure  process  analysis  of  a  real  system 


In  order  to  verify  that  the  assumptions  stated  in  Section  2  lead  to  a  better  modeling  of  failure 
processes  than  other  models,  an  experiment  was  designed.  The  experiment  consisted  in  the 
adquisition  of  data  concerning  both  the  failure  process  and  the  use  of  a  general  purpose  time  sharing 
system.  The  system  choosen  was  the  CMU-A,  a  PDP-10  used  by  the  Computer  Science  Department  at 
Camegie-Mellon  University  as  the  main  general  purpose  computational  tool.  The  system  consists  of  a 
KL-10  processor,  one  megaword  of  memory,  eight  disk  drives  totalling  1600  megabytes  of  online 
storage  and  two  magnetic  tape  drives.  The  system  runs  a  slightly  modified  version  of  the  standard 
TOPS- 10  operating  system  [Bell  78aJ. 

The  software  packages  used  to  instrument  the  experiment  are  illustrated  in  Figure  3-1 .  Information 
about  failures  is  obtained  from  an  online  error  log  file  maintained  by  a  system  program,  which  records 
the  information  produced  by  different  error  formatting  routines.  Entries  are  made  to  this  file  for  each 
hardware  error  detected  in  the  system,  for  system  reloads,  for  disks  performance  statistics,  and  so  on 
[Digital  78].  The  error  log  is  later  processed  by  SEADS,  a  FORTRAN  package  which  allows  to  list  the 
times  of  detection  of  errors  associated  with  a  particular  resource.  In  order  to  obtain  accurate 
information  about  the  use  of  the  system,  a  special  SAIL  program,  SYSMON,  was  written  that  samples 
the  values  of  30  system  parameters  twice  every  five  minutes,  the  two  samples  in  a  five  minutes  interval 
being  one  second  apart  In  this  way,  I/O  traffic,  system  overhead  values,  etc.,  can  be  obtained 
averaged  on  a  one  second  interval  or  in  a  five  minute  interval  with  a  resolution  of  5  minutes.  The  files 
generated  by  SYSMON  are  later  processed  by  another  SAIL  package,  READSY,  which  computes  the 
periodic  component  and  autocorrelation  function  of  the  utilization  function  of  a  particular  system 
resource.  The  information  generated  by  SEADS  and  READSY  is  then  processed  by  an  APL  package 
(POWELL)  which  estimates  the  maximum  likelihood  parameters  of  the  pdf  of  the  time  to  failure  of  a 
particular  resource.  Finally,  in  a  separate  SAIL  package,  C2TST,  the  values  predicted  by  the 
cyclostationary  model  and  other  models  described  in  Section  4  are  compared  with  the  information 
stored  in  the  error  log  according  to  a  x2  goodness-of-fit  test 

The  operational  policies  regulating  the  use  of  this  system  at  CMU  make  it  a  good  starting  point  to 
check  the  validity  of  the  ideas  exposed  in  Section  2.  Its  steady  state  operation  during  weekdays  can 
^  be  understood  from  Figures  2-4  and  2-7.  Recall  that  this  figure  plots  the  sampled  values  of  the  fraction 
V  of  time  considered  to  be  operating  system  overhead  for  five  consecutive  weekdays.  The  value  of  the 
accumulated  overhead  time  is  obtained  by  executing  a  Monitor  Call  and  includes  the  time  spent  in 
clock  queue  processing,  short  command  processing,  swapping  and  scheduling  decision,  and 
software  context  switching  [Digital  77].  This  value  does  not  include  Monitor  Calls  execution  nor  I/O 
interrupt  times.  It  is  not  exactly  the  time  that  the  system  is  executing  in  kernel  mode,  but  it  is  close 
enough  for  our  purposes. 
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Figure  3-1:  Software  packages  used  in  the  validation  of  the 
cyclostationary  modeling  methodology. 
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3.1  Probability  Distribution  Function  of  the  Time  to  Failure  of 
a  File  System 

Figure  2-2  shows  the  results  of  compiling  five  days  of  disk  utilization  samples  into  a  single  24  hour 
period.  Along  with  the  estimated  average,  this  figure  shows  the  function  mdk(t)  obtained  from  a  finite 
Fourier  series  expansion  (see  Appendix  I  for  details).  A  Fourier  series  expansion  is  a  least  squares  fit 
to  our  data  and  is  a  good  way  of  eliminating  the  "noise”  present  in  the  estimated  average  due  to  the 
finiteness  of  the  sample.  The  data  in  Figure  2-2  corresponds  in  fact  to  the  function  mdk(t)  in  Section 
2.3.  after  sampling  its  values  every  five  minutes.  After  substracting  from  b(t)  the  value  of  mdk(t),  the 
sampled  values  of  the  process  zdk(t)  are  available  for  estimation  of  its  autocorrelation  function. 


Figure  3-2:  Estimated  and  approximated  Autocorrelation 
functions  of  the  file  system  utilization  process. 

Figure  3-2  shows  the  estimated  autocorrelation  function  for  the  process  zdk(t).  From  its 
appearance,  it  seems  that  an  autocorrelation  function  of  the  form 


Ra(t) 


«ie 


a2e 


(3.1) 


would  be  appropiate  to  approximate  the  real  autocorrelation  function.  The  noisy  appearance  of  the 
estimated  autocorrelation  function  is  again  a  consequence  of  the  finite  sample  size  available  for  its 
computation.  The  main  problem  in  the  evaluation  of  the  a,  and  /J.  is  that  they  are,  in  principle,  very 
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sensitive  to  the  sampling  interval  (in  this  case,  5  minutes).  In  the  Appendix  III,  the  exact  procedure 
followed  to  evaluate  them  is  described  in  detail. 

With  the  autocorrelation  function  given  in  (3.1),  the  following  expression  is  obtained  for  the 
variance  o^t) : 


<r2( t)  »  2a,  f  (t-r)e^1  dr  +  2a2  f  (t-r)e^2  dr 


‘  T  2a, 


^2J  j82 

and  substituting  (3.3)  in  (2.23)  we  obtain, 


n 


„  „  W  CTdk1 .  ‘^1*.  gdk2  .  ’^1 

"  V  ‘  h  1  1 

where  the  following  constants  have  been  defined 


“dk  *  ®dk*^dk  +  Cdk  '  ^dk^Sdk,bmin^ 


ffdk1  *  ■^"Sdk 

gdk2  »  ^4 


(3.2) 


(3.3) 


(3.4) 


(3.5) 

(3.6) 

(3.7) 


The  hazard  function  is  given  by 


hdk<T>  *  «dk 


i  a*dk<T> 
^dk(T) 


(3.8) 


The  statistics  of  the  time  to  failure  for  a  doubly  stochastic  Poisson  process  when  the  intensity 
process  is  a  cyclostationary  process  are  then  equivalent  to  the  statistics  of  a  non  homogeneous 
Poisson  process  with  hazard  function  given  in  (3.8).  Although  impressive,  this  hazard  function 
reduces  to  a  constant  term  plus  a  periodic  component  plus  an  exponentially  decreasing  term.  Note 
that  neglecting  the  periodic  component,  this  hazard  function  is  exponentially  decreasing  with  the 
following  extreme  values 
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Figure  3*3:  Hazard  function  of  the  equivalent  non  homogeneous 
Poisson  process  characterizing  the  statistics  of  the  time  to  failure 
of  a  file  system. 


(3.9) 
(3.10) 

as  shown  in  Figure  3-3. 


hdk(0)  *  «dk 

hdk(»)  =  adk-adk1-<rdk2 


3.2  The  Probability  Distribution  Functions  of  the  Time  to 
System  Failure 


The  periodic  component  of  the  kernel  utilization  process,  m/t),  has  been  shown  in  Figure  2*5. 

*7 

Figure  3-4  shows  the  autocorrelation  function  of  the  process  z^t),  suggesting  again  an 
approximation  of  the  form  given  in  (3.1 ).  The  following  expression  is  then  obtained  for  the  PDF  of  the 
time  to  system  failure 


P^Kr) 


4a  . 
l-^(r)e(  *  * 


(3.11) 


where 


3  6  9  12  IS  18  21  24  27  30  33  36  39  42  45  48 

Time  (hours) 


Figure  3*4:  Estimated  and  approximated  Autocorrelation 
functions  of  the  kernel  utilization  process. 


®sy  *  +  ®hw^sy  +  ®hw  '  Psy(®sw  +  ®hw’^min^  *  ^ay^aw'N)^ 

®sy1  “  ^Ssw  +  3hw^ 

®sy2  “  p2  +  ®hw^ 


(3.12) 

(3.13) 


(3.14) 


The  hazard  function  is  given  by 


VT)  "  <VffsyX1  -Jl  -e^] 


1  a^(r) 

^(T)  dr 


(3.15) 


3.3  Simplified  expressions  for  known  starting  time 


All  the  expressions  given  in  Sections  3.1.  and  3.2.  have  been  obtained  after  computing  the 
expectation  for  ail  possible  values  of  the  starting  time  in  a  one  day  period.  If  the  system  starting  time  is 
known,  different  expressions  are  obtained.  The  only  differences  between  the  PDF  of  the  time  to  failure 
with  known  starting  time  and  the  PDF  averaged  over  a  one  day  period  is  that  the  function  4(r) 
becomes  a  constant  equal  to  one  and  that  the  a  term  in  the  exponential  is  slightly  different.  In 
particular,  for  the  case  of  the  file  system  failures. 
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P-vtKrjts) 


, .  e  “'dk  ■  <«V»dM-W  • 


(3.16) 


where  Pdk(Kr|ts)  is  the  probability  that  a  failure  will  be  detected  before  time  r  +  ts  given  that  no  failure 


had  bean  detected  at  time  ts,  and 


®  dk  *  s<jk^dk^T  +  ^  (3.17) 

“  ”dk  *  Cdk  *  Pdk^dk'^rrtJ  (3.18) 

*  f  Ridk(ir)6T  (3.19) 

J0 

The  hazard  function  for  known  starting  time  is 

bdk(Tlts)  *  Sdkmdk^*s  +  T^  +  “  ’dk'^dklt1  "e  1  3  '  “awO  "®  (3.20) 


Similar  expressions  can  be  derived  for  the  distribution  of  the  time  to  system  failure. 


3.4  Distribution  functions  of  the  time  to  system  failure  due  to 
software  and  of  the  time  to  system  failure  due  to 
hardware  transients 


Once  the  values  of  *w,  s^,  c^,  kg  are  known,  it  is  straightforward  to  derive  an  expression  for  the 
PDF  of  the  time  to  system  failure  due  to  hardware  transients.  Repeating  the  derivation  described  in 
Section  2.4.  with  sw  »  kQ  »  0,  the  following  expression  is  obtained 


PkJKt) 


. . *JT, 


(3.21) 


where 
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ahw  *  ®hw^sy  +  Chw  *  Psy^hw^min)  (3*22) 

ahw1  *  ^"*h w  <3-23) 

ffhw2  *  jj*  ®hw  (3.24) 

In  the  general  case,  to  obtain  the  PDF  of  the  time  to  system  failure  due  to  software  errors,  a  similar 
equation  to  (3.16)  would  be  obtained,  but  with  the  following  parameters 


asw  *  Ssw^*sy  *  Pay^sw’^min^  *  Fsy^sw'^O^ 

a  ,  -  -^s2 

swl  a  sw 
®»  0 

CTsw2  *  ®w* 


(3.25) 

(3.26) 

(3.27) 


Note  that  if  the  system  failure  processes  due  to  software  errors  and  hardware  transients  are 
considered  to  be  nonhomogeneous  Poisson  processes,  each  with  a  PDF  of  the  form  (3.16),  the 
superposition  of  both  processes  (i.e.,  the  process  obtained  by  adding  the  hazard  functions  of  the 
software  failure  process  and  the  hardware  failure  process)  is  not  equal  to  the  total  system  failure 
process,  whose  PDF  is  given  by  (3.11).  This  is  because  they  are  not  statistically  independent.  Indeed, 
both  failure  processes  have  a  common  cause,  the  utilization  process  of  the  kernel  of  the  operating 
system. 


Table  3-1  gives  the  maximum  likelihood  values  of  sM,  s^,  c^,  and  kg  for  the  CMU-A,  along  with  the 
value  of  Note  that  since  the  value  of  Kg  is  larger  than  (Tt^,  expressions  (3.22)  thru  (3.25)  may  not 
be  valid.  The  correction  term  p(saw,kQ)  has  been  computed  assuming  that  ^<<<1*^,  condition  that  the 
maximum  likelihood  value  of  Kg  doss  not  verify.  In  fact,  if  KgXh^  the  probability  density  function  of 
the  time  to  failure  due  to  a  software  error  degenerates  into  an  exponential,  such  that 

p«*(Kt>  ■  1  -  e  P*y<**^°)T  (3.28) 


the  PDF  of  the  time  to  system  failure  due  to  hardware  transients  being  given  in  (3.21). 

Figure  3-5  shows  the  relationship  between  the  instantaneous  value  of  the  system  failure  rate  and 
the  software  and  hardware  components.  Note  how  the  software  failure  rate  is  zero  for  a  wide  range  of 
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Parameter 

Value 

8«w 

0.158 

®hw 

0.086 

^hw 

0.0079 

c 

0.225 

iTl 

0.19 

«y 

Table  3-1:  Maximum  likelihood  values  of  the  coefficients 
defining  the  relationship  between  kernel  utilization  and  system 
failure  rate. 

values  of  k(t),  but  that  its  slope  is  larger  than  the  slope  of  the  failure  rate  due  to  hardware  errors. 
Figure  3-5  thus  suggests  that  to  assume  a  linear  relationship  between  the  system  failure  rate  due  to 
software  errors  and  kernel  utilization  may  be  an  oversimplification.  In  fact,  it  seems  reasonable  to 
expect  the  probability  of  observing  a  software  error  to  increase  with  the  length  of  time  that  the 
software  is  exercised  and  with  a  "stress”  factor  depending  on  the  apparent  complexity  of  the  input 
data  to  be  processed  at  a  given  time.  In  this  case,  perhaps  a  higher  degree  polynomial  would  better 
describe  the  relationship  between  the  software  failure  rate  and  software  utilization. 


k(t) 

Figure  3-5:  Relationship  between  the  system  failure  rate  and  its 
software  and  hardware  components. 
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4.  Discussion 


Although  the  assumptions  made  in  Section  2  are  less  restrictve  than  the  usual  assumption  of 
modelling  the  failure  process  with  a  constant  failure  rate,  the  validity  of  the  methodology  presented 
here  can  be  asserted  only  by  comparison  with  the  behavior  of  a  real  system  and  contrasting  the 
results  that  would  be  predicted  by  traditional  models.  This  is  the  subject  of  Section  4.1.,  where  the 
results  given  in  Section  3  are  compared  with  the  values  predicted  by  assuming  either  an  exponential 
distribution,  a  Weibull  distribution,  and  a  periodic  distribution  for  the  time  to  failure.  In  section  4.2.  an 
explanation  is  given  for  why  the  apparent  failure  rate  is  decreasing,  and  finally  in  Section  4.3.  some 
preliminary  conclusions  are  summarized. 

4.1  Comparisons  with  other  models 

The  more  widespread  model  used  to  characterize  the  failure  process  of  digital  computers  assumes 
the  failure  process  to  be  a  homogeneous  Poisson  process.  The  PDF  of  the  time  to  failure  is  then  given 
by 


Pe(Kr) 


1  -  e 


(4.1) 


where  Xe  is  the  (constant)  failure  rate.  The  maximum  likelihood  estimate  of  \s  is  obtained  simpy  by 
dividing  the  time  that  the  system  has  been  operational  by  the  number  of  failures  reported.  All 
functions  and  parameters  related  to  this  model  will  be  noted  with  subindex  He"  and  from  now  on  this 
model  will  be  referred  to  as  the  exponential  model. 

However,  empirical  studies  [McConne!  79a],  [Wagoner  73]  have  shown  that  a  Weibull  distribution 
gives  much  better  goodness  of  fit  to  experimental  data  than  a  simple  exponential.  The  Weibull  PDF  is 
given  by 


Pw(Kt)  «  1-e 


(Xwr)aw 


(4.2) 


The  Weibull  distribution  is  characterized  by  two  parameters :  Xw  ,  the  scale  parameter,  and  aw  ,  the 
shape  parameter.  For  aw  ■  1 ,  the  Weibull  distribution  degenerates  to  the  exponential.  For  <*W>1 ,  the 
Weibull  distribution  has  an  increasing  failure  rate.  A  decreasing  failure  rate  corresponds  to  V1A" 
reports  published  to  date  claim  that  a  decreasing  failure  rate  Weibull  distribution  fits  experimental 


40 


WORKLOAD,  PERFORMANCE,  AND  RELIABILITY  OF  DIGITAL  COMPUTING  SYSTEMS 


data  much  better  than  a  plain  exponential  model.  Numerical  procedures  have  been  developed  to  find 
the  maximum  likelihood  estimates  of  Aw  and  aw.  These  procedures  are  based  on  the  works  of 
[Thoman  69,  Berger  74,  Romano  77]  and  FORTRAN  programs  implementing  them  are  given  in 
[McConne!  79b]. 

A  workload  dependent  model  has  been  presented  in  [Butner  80].  A  linear  dependency  between 
failure  rate  and  workload  is  also  assumed.  The  workload  is  characterized  by  a  periodic  function  of 
time.  The  PDF  becomes  an  exponential H modulated”  by  a  periodic  function 

-K  t  -f  u  (T) 

Pp(Kr)  -  1  -e  p  e  p  p  (4.3) 

where  Fp  is  defined  as  the  load  induced  failure  rate  and  U(r)  denotes  the  instantaneous  load  value. 
This  model  will  be  referred  to  as  the  periodic  model,  ail  its  parameters  having  the  subindex  Hp".  Using 
the  notation  developed  in  Section  2,  this  is  equivalent  to  assume  an  utilization  function  u(t)  ■  m(t), 
where  only  the  periodic  component  is  taken  into  account,  and  where  the  Gaussian  process  z(t)  has 
been  neglected.  In  this  case, 

p(t)  -  E{\p(t)e  'j£  Xp(T)dT  }  (4.4) 


where 

V*)  -  spm(t)  +  cp 


(4.5) 


and 


P(Kr) 


•<» m  *  c  Jt  In  #T) 

1  -e  p  *  e 


(4.6) 


Note  that  (4.3)  and  (4.6)  are  equivalent.  In  Section  11.5.  the  equations  for  computing  the  maximum 
lykelihood  values  of  sp  and  cp  from  a  history  of  failures  are  given. 

Finally,  the  model  presented  in  this  report  will  be  referred  as  the  cyclostationary  model.  An 
expression  for  its  PDF  is  rewritten  here 
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Figure  4-1:  Two  alternatives  to  characterize  system  reliabiHty. 
The  maximum  likelihood  values  of  the  hazard  function 
parameters  can  be  evaluated  from  the  resource  utilization 
functions  or  directly  from  a  history  of  failures. 
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Table  4.1.  summarizes  the  densities,  Reliability  and  Hazard  functions  for  each  of  the  four  models. 
Note  that  the  cyclostationary  model  has  both  an  asymptotically  decreasing  failure  rate  and  a  periodic 
component.  Qualitatively,  the  cyclostationary  model  seems  to  integrate  the  approaches  of  the  Weibull 
and  periodic  models.  Note  also  that,  for  the  case  of  a  file  system  failures,  only  two  parameters  need  to 
be  estimated  from  a  history  of  system  failures  (sdk  and  cdk),  the  other  parameters  in  equation  (4.7) 
being  measured  from  the  actual  system  behavior  (the  periodic  component  and  the  autocorrelation 
function  of  the  resource  utilization  process).  Since  the  cyclostationary  model  suggests  a  PDF  of  the 
form  shown  in  (4.7),  it  is  conceivable  to  postulate  (4.7)  as  the  real  PDF  of  the  failure  process  and 
estimate  the  values  of  ac,  ac1,  <rc2,  fiv  p2  directly  from  a  history  of  failures,  therefore  avoiding  the 
measurement  of  the  resource  utilization  functions.  Figure  4-1  describes  these  two  alternatives 
available  to  characterize  system  reliability.  In  Section  11.3.  the  equations  used  to  estimate  the  values 
of  these  parameters  drirectly  from  a  history  of  system  failures  are  given. 

The  fifth  distribution  in  Table  4.1.  is  a  simplified  version  of  the  distribution  obtained  with  the 
cyciostationary  model,  considering  only  one  exponential  in  the  hazard  function,  and  neglecting  the 
periodic  component  <P(t).  Section  11.4.  gives  the  equations  for  estimating  the  maximum  likelihood 
parameters  of  this  distribution  from  a  history  of  system  failures. 

Next,  quantitative  comparisons  using  data  of  a  real  system  are  in  order.  Table  4-2  show  the  results 
of  applying  a  x2  goodness  of  fit  test  between  the  actual  failures  observed  in  the  file  system  described 
in  section  3.1 .  and  the  distributions  predicted  by  the  above  four  models.  A  x2  value  smaller  than  x^ 
(i.e.,  a  level  of  confidence  greater  than  0.05)  indicates  a  good  fit  between  predicted  and  observed 
behavior  and  suggests  the  acceptance  of  the  hypothetical  distribution  as  the  real  distribution 
underlying  the  failure  process. 

As  can  be  seen  from  Table  4-2,  only  the  cyciostationary  model  (both  with  direct  and  indirect 
evaluation  of  its  maximum  likelihood  parameters)  show  a  good  fit  with  experimental  data.  Neither  the 
exponential  nor  the  periodic  failure  rate  models  seem  to  be  able  to  describe  the  failure  process  with 
significant  accuracy.  The  simplified  cyciostationary  model  distribution  and  the  Weibull  distibution  are 
almost  in  the  border  line  of  acceptance.  Further  insight  can  be  gained  by  direct  comparison  of  the 
hazard  functions  of  the  above  four  models.  Figure  4-2  shows  the  hazard  function  of  the  above  four 
models  for  the  case  of  file  system  failures. 
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Table  4-3  shows  the  results  of  applying  a  x2  goodness-of-fit  test  to  the  four  models  in  the  case  of 
system  failures.  Note  that  for  the  periodic  model,  the  maximum  likelihood  values  of  the  coefficient  are 
such  that  the  proportionality  term  vanishes,  and  the  constant  term  equals  the  \  value  of  the 
exponential  model.  Although  all  models  give  a  level  of  significance  larger  that  0.05,  the 
cyclostationary  model  is  again  clearly  superior  giving  levels  of  significance  of  0.9.  The  hazard 
functions  of  the  four  models  for  the  case  of  system  failures  are  given  in  Figure  4-3. 

4.2  The  decreasing  hazard  function  paradox 

Tha  hazard  function  found  in  the  cyclostationary  model  presents  the  following  paradox  :  neglecting 
the  periodic  component,  expression  (3.15)  means  that  no  matter  at  what  time  we  start  observing  a 
system,  the  statistics  of  the  time  to  failure  are  equivalent  to  the  statistics  of  a  non  homogeneous 
Poisson  process  with  decreasing  hazard  function.  Since  the  hazard  function  is  roughly  the  rate  at 
which  failures  will  be  detected,  this  means  that  no  matter  at  which  time  we  start  observing  a  system 
the  apparent  rate  at  which  failures  are  detected  will  be  a  decreasing  funtion  of  time.  In  this  section, 
the  reason  for  such  surprising  behavior  is  investigated  and  explained. 

To  understand  the  decreasing  hazard  function  paradox,  start  with  the  simplest  possible  case. 
Assume  that  the  real  failure  rate  is  given  by  a  constant  plus  white  noise. 


Aa1(t)  *  m  +  x,(t) 


(4.18) 


where  m  is  the  (constant)  mean  failure  rate  and  x.,(t)  is  a  stationary,  zero  mean  Gaussian  process  with 
autocorrelation  function 


w. 


Rx  x  (r)  -  -^<r) 
x1*1  2 


(4.19) 


and 


6(t) 


Rm 

h<l>0 


1/h  0<T<Jh 
0  r>h 


Assume  that  m»Wt  such  that  the  probability  of  X#1(t)  being  negative  can  be  neglected.  The 
probability  density  function  of  the  time  to  failure  is  then  given  by 
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exponential 


R.(t)  -  e 

h#(T)  -  X. 


-X.T 


WeibuW 


Rw(r)  *  0 


(X^Tj^W 


hw(r) 


(Xwt) 


la 


(4.8) 

(4.9) 


(4.10) 

(4.11) 


Periodic 


Rp(r)  -  .X>V'«n 
*  [\  *  ] 


Cvclostationarv 


Rc(r)  -  e 


-fl.T  O  -  -/J,T 

-<Xc  +  ffc1  +ac2>T-^LI1e  1  l1'*  2  I  + 


Simplified  Cvclostatlonafv 

Y  -fl  r 
-<amYIJT-*sa(i  «  m  1 
Rm(r)  *  e  mlm>  Pm 


(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 


Table  4*1 :  Reliability  and  Hazard  functions  of  the  five  compared 
models. 


DISCUSSION 


46 


Model 

Parameter 

Values 

Degrees  of 
Freedom 

X2  value 

*0.05 

Level  of 
Confidence 

Exponential 

\e  -0.67 

7 

130 

14.067 

0 

Weibull 

Aw  -0.91 
*0.68 

8 

17.717 

15.507 

0.026 

Periodic 

s  -1.25 

Cp  =0.28 

12 

1007.194 

21.026 

0 

Cyclostat. 

Sc  -  14.00 

Cc  -  2.01 

8 

8.69 

15.07 

0.36 

Cyclostat. 

(Direct) 

ac  -2.13 

ffic  =1-42 

a2c  =  403 
)51  -0.59 

02  -0.21 

6 

8.642 

12.592 

0.19 

Simplified 

Cyclostat. 

ac  -1.69 
a^-1.38 
/91  =1.38 

8 

19.434 

15.507 

0.013 

Table  4-2:  Results  of  a  x2  goodness-of-fit  test  with  the 
Exponential,  Weibull,  Periodic,  and  Cyclostationary  models  for  file 
system  failures.  Only  the  Cyclostationary  model  gives  a  level  of 
confidence  greater  than  0.05  The  Weibull  and  simplified 
cyclostationary  models  give  smaller  levels  of  confidence  but  close 
to  0.05.  The  hypothesis  that  the  time  to  failure  can  be 
characterized  with  Exponential  or  Periodic  models  has  to  be 
rejected.  The  data  used  was  obtained  from  five  weekdays  of 
system  operation  during  which  877  (transient)  failures  were 
detected.  The  MTTF  value  is  7  minutes.  The  file  system  is 
composed  of  8  RP06  disk  drives  totalling  1600  megabytes  of  on 
line  storage. 
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Figure  4*2:  Hazard  functions  predicted  by  Exponential,  Weibull, 
Periodic,  and  Cyclostationary  models  for  file  system  failures. 
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Model 


Parameter 

Values 


Degrees  of 
Freedom 


X2  value 


*0.05 


Level  of 
Confidence 


Exponential 


\a  *  0.0073 


7.87 


15.507 


0.45 


Weibull 


\w  *  0.0074 
aw  *0.98 


7.95 


14.067 


0.35 


Periodic 


s  =0.0 
Cp  =  0.0073 


Cyclostat. 


Ssw  a  0.158 
shw  =  0.0869 
Chu(  =  0.0079 
k0  a  0.0357 


1.61 


11.070 


0.9 


Cyclostat. 

(Direct) 


a  a  0.013 
<r1c  »  0.0054 
a*.  =0.0080 
*  0.21 

0,  a  0.0041 


1.66 


11.070 


0.9 


Simplified 

Cyclostat. 


a  =0.014 
a.  =0.0064 
0,  =0.21 


0.75 


12.592 


0.9 


Table  4*3:  Results  of  a  x2  goodness-of-fit  test  with  the 
Exponential,  Weibull,  Periodic,  and  Cyclostationary  models  for 
system  failures  (crashes).  Although  all  models  give  a  level  of 
confidence  larger  than  0.05,  the  Cyclostationary  model  shows  a. 
better  fit  to  real  data.  Note  that  for  the  Periodic  model  the 
maximum  likelihood  values  of  the  coefficients  is  such  that  the 
proportionality  coefficient  vanishes  and  the  constant  term  is  equal 
to  the  X  value  of  the  Exponential  model.  The  data  used  was 
obtained  from  29  weekdays  of  system  operation  during  which  60 
failures  were  detected  giving  a  MTTS  (Mean  Time  To  restart) 
value  of  1 1  hours. 
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Figure  4*3:  Hazard  functions  predicted  by  Exponential,  Weibull, 
Periodic,  and  Cyclostationary  models  for  system  failures. 
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where  <r2(t)  is  the  variance  of  the  integrated  process 

<r*(t)  -  W,  /  (t-r)  5(r)  dr 

/o 

■  W1t 


such  that 


>.«  ■ 


•rot  W.t  \ 

e  e  1  ; 


P(Kt)  -  l-e 


-<m>w  1)r 


hn(r)  »  m-W1 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


(4.24) 


(4.25) 


This  failure  process  is  equivalent  to  a  homogeneous  Poisson  process  with  an  apparent  hazard 
function  not  equal  to  the  mean  failure  rate,  but  equal  to  the  difference  of  the  mean  failure  rate  minus 
the  "power"  of  the  noise,  Wr  The  reason  for  that  is  that  the  failure  rate  appearing  in  the  exponent  of 
an  exponential,  variations  above  the  mean  failure  rate  are  not  equally  weighted  with  variations  below 
the  mean.  In  fact,  the  variations  below  the  mean  value  are  more  heavily  weighted,  and  hence  the 
resulting  smaller  limiting  failure  rate  when  the  expectation  is  taken  over  all  possible  realizations  of  the 
failure  rate  process. 

Assume  now  that  the  real  failure  rate  is  equal  to  a  constant  plus  a  zero  mean,  stationary  Gaussian 
process 


X2(t)  •  m  +  xa(t) 


(4.26) 


but  that  now  the  autocorrelation  function  of  the  Gaussian  process  is  given  by 
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(T) 


JL 

2 


-0M 

e 


In  this  case, 

Defining 


w  3  JL 
2  0 


the  following  expressions  are  obtained 

w,  -0t 

a  /  -mt  W0t  •  — 1  * «  ]\ 

P2(t)--^-(e  e  2  0  J) 

w,  fir 
-{m-WJr  •  -~f1  •  e  J 
P2(Kr)-1-e  *  0 


h2(r)  »  m-W2[l  -e‘^T] 


(4.27) 


(4.28) 


(4.29) 


(4.30) 


(4.31) 


(4.32) 


This  failure  process  is  then  equivalent  to  a  non  homogeneous  Poisson  process  with  an 
exponentially  decreasing  hazard  function.  For  r«0,  the  apparent  hazard  function  is  equal  to  the 
mean  real  failure  rate,  m.  For  r  ■  oo,  the  apparent  hazard  function  equals  the  same  value  that  had 
been  obtained  assuming  the  Gausian  process  to  be  white  noise.  And  as  r  increases,  the  failure  rate 
approaches  this  limiting  value  exponentially. 

Finally,  note  that  if  W1  -  W2,  the  system  with  white  noise  utilization  process  will  be  more  reliable 
than  another  system  having  a  utilization  process  with  autocorrelation  function  given  by  (4.27).  In  the 
case  of  white  noise,  the  system  reaches  the  minimum  value  of  its  hazard  function  at  t  ■  0,  while  for  an 
autocorrelation  function  of  the  form  (4.27)  the  minimum  value  is  approached  exponentially.  Hence  a 
non  obvious  way  of  increasing  the  reliability  of  a  particular  resource  would  be  to  build  a  system  such 
that  the  utilization  process  of  that  resource  approaches  white  noise  as  much  as  possible. 
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4.3  Preliminary  conclusions 

It  has  been  shown  how  the  cyclostationary  model  is  capable  of  predicting  the  reliability  of  a  Time¬ 
sharing  system  in  steady  state  operation.  Workload  dependency  is  stated  explicitely  in  the  model,  by 
means  of  resource  utilization  functions.  System  reliability  is  evaluated  in  terms  of  utilization  of  system 
singularities  (i.e.,  the  Kernel  of  the  operating  system).  And,  in  general,  resource  reliability  is  evaluated 
in  function  of  the  utilization  of  each  resource. 

Sofware  and  hardware  reliability  can  be  evaluated  separately  merely  by  observing  a  history  of 
system  failures  and  some  knowledge  of  how  the  system  behaves  (periodic  mean  and  autocorrelation 
function).  A  linear  relationship  between  software  failure  rate  and  software  utilization  has  given 
somewhat  contradictory  results,  suggesting  that  perhaps  more  complex  relationships  need  to  be 
considered.  In  any  case,  it  has  been  shown  how  establishing  the  relationships  between  software 
failure  rate,  hardware  failure  rate,  and  kernel  utilization,  it  is  in  principle  possible  to  evaluate  the 
contribution  of  software  and  hardware  to  the  unreliability  of  the  total  system. 

Perhaps  one  of  the  more  important  results  is  that  the  probability  density  function  for  the  simplified 
cyclostationary  model  (having  a  single  exponential  in  its  hazard  function)  has  a  known  Laplace 
transform,  making  it  suitable  for  Markov  modelling.  Neither  the  complete  Cyclostationary,  nor  the 
Weibuil,  nor  the  Periodic  models  lead  to  probability  density  functions  with  known  Laplace  transforms. 

From  a  user  viewpoint,  there  is  a  reinforcement  effect  between  workload  and  lack  of  reliability. 
Higher  workload  implies  that  the  Kernel  of  the  operating  system  has  to  take  more  decisions  per  unit 
time,  increasing  the  probability  of  a  system  failure.  Hence,  not  only  the  user  receives  less  CPU  cycles 
per  unit  time,  but  the  probability  that  these  cycles  will  become  useless  because  the  job  will  have  to  be 
restarted  also  increases. 

Hence,  high  reliability  seems  to  be  in  contradiction  with  other  performance  measures  (such  as  the 
maximum  number  of  jobs  allowed  to  be  simultaneoulsy  active). 

But  the  contradiction  between  reliability  and  other  performance  measures  seems  to  be  of  a  deeper 
nature.  In  [Spim  77]  several  paging  algorithms  are  described  and  modelled.  Page  faulting  in  a  virtual 
memory  system  can  be  described  as  a  stochastic  process,  and  a  usual  optimality  criteria  for  paging 
algorithms  is  how  well  they  are  able  to  perdict  future  system  behavior  given  past  and  present  system 
behavior.  This  is  exactly  the  information  given  by  the  autocorrelation  function,  and  in  [Spim  77] 
several  paging  algorithms  are  compared  in  terms  of  how  well  their  predictions  fit  the  autocorrelation 
function  of  a  real  page  faulting  process.  But  in  Section  4.2.  it  has  been  shown  that  a  way  of  improving 
reliability  is  to  have  white  noise  as  the  resource  utilization  process.  For  white  noise,  the  future  values 
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of  the  process  are  completelly  unpredictable  no  matter  for  how  long  has  the  system  been  observed 
and  no  matter  how  close  in  the  future  is  the  prediction  desired.  Hence,  for  an  optimally  reliable  paging 
system,  its  utilization  process  should  be  white  noise.  Since  future  system  behavior  would  be 
unpredictable,  an  optimum  paging  algorithm  under  these  conditions  would  just  swap  out  of  memory 
pages  at  random,  therefore  lowering  system  performance. 
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5.  Proposed  research 


5.1  On  the  linear  dependency  between  overhead  and  failure 
rate 

The  first  topic  to  be  investigated  in  depth  would  be  the  real  dependency  between  overhead  and 
failure  rate.  A  linear  relationship  has  been  assumed  for  the  cyclostationary  model,  and  any  other 
relationship  may  lead  to  hopelessly  complex  mathematical  problems  in  the  evaluation  of  the 
expectation  of  the  modified  Gaussian  process.  However,  there  is  always  a  possibility  that  other 
dependencies  may  be  more  accurate,  and  even  if  exact  expressions  for  the  distribution  of  the  time  to 
failure  cannot  be  obtained  for  dependencies  other  than  linear,  errors  due  to  a  linear  relationship 
assumption  should  be  understood. 

Since  a  failure  cannot  be  detected  if  the  system  is  not  used,  the  only  a  priori  assumption  that  seems 
reasonable  is  that  the  failure  rate  must  be  a  non-decreasing  function  of  the  system  overhead.  What 
needs  to  be  known  is  for  what  ranges  a  linear  dependency  is  accurate,  what  are  the  confidence 
intervals  that  can  be  obtained,  and  to  explore  the  possibility  of  characterizing  the  failure  rate  with 
relationships  other  than  linear  if  necessary. 


5.2  Generalization  to  systems  showing  a  non-cyclostationary 
behavior 

One  of  the  fundamental  assumptions  made  to  develop  the  cyclostationary  model  has  been  that  the 
system  overhead  could  be  approximated  by  adding  a  periodic  function  to  a  modified  Gausian 
process.  This  may  be  a  good  approximation  for  time-sharing  systems,  but  it  certainly  does  not  apply  to 
many  real-time  and  command  and  control  systems.  In  fact,  the  highest  demand  for  high  availability 
systems  comes  from  special  purpose  command  and  control  systems  like  the  ones  to  be  installed  in 
aircrafts,  missiles,  satelites,  and  so  on.  For  some  of  these  systems,  the  workload  can  be  modelled  by  a 
sequence  of  load  states.  If  the  exact  sequence  is  not  known  in  advance  but  the  possible  alternatives 
are  known,  the  instantaneous  value  of  the  mean  workload  could  be  modeled  by  a  aemi-markov 
process  [Howard  71], 

The  cyclostationary  model  would  then  evolve  to  a  model  in  which  the  instantaneous  mean  failure 
rate  is  not  a  periodic  function  of  time,  but  a  random  variable  whose  statistics  depend  on  the  mission  to 
accomplish.  This  new  model  would  have,  in  addition,  the  ability  to  incorporate  the  effect  of  permanent 
hardware  failures,  transient  hardware  failures,  and  software  failures.  In  fact,  most  performance- 
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reliability  models  presented  to  date  just  assume  that  in  the  presence  of  a  permanent  hardware  failure 
the  system  reconfigures  itself  and  continues  to  operate  with  a  possibly  different  computational 
capacity.  Whether  the  capacity  diminishes  or  the  workload  increases,  the  result  will  be  a  system 
operating  in  one  of  a  set  of  possible  states  for  a  given  period  of  time. 


5.3  Characterizing  total  system  performance 

Present  reliability  evaluation  tools,  such  as  Reliability  or  Availability,  are  felt  to  be  innadaquate  due 
to  the  large  gap  that  separates  say,  the  Availability  of  a  computing  facility  and  the  cost  that  has  to  be 
paid  due  to  of  lack  of  reliability.  Digital  computers  are  used  to  store  and  process  information,  and  the 
sooner  the  desired  information  is  available,  the  better  the  system.  The  occurrance  of  a  system  failure 
means  waiting  until  operation  is  restored,  bringing  the  machine  to  a  consistent  state,  possibly 
restarting  computations  that  were  interrupted  because  of  the  failure  and  (if  possible)  updating  the 
system  with  the  information  it  was  supposed  to  process  while  it  was  not  operational.  In  short,  it  means 
a  delay  in  obtaining  the  desired  information  and  an  added  cost  associated  with  the  extra 
computations  related  to  restoring  the  system  to  the  desired  state  after  the  failure  occurred. 

From  a  single  user  viewpoint,  a  failure  also  means  a  delay.  In  [Castillo  80]  the  expected  eilapsed 
time  required  to  complete  a  program  was  computed  under  rather  restrictive  assumptions,  but 
separating  the  "useful"  time  that  leads  to  program  completion  from  the  "useless"  time  due  to  lack  of 
reliability.  Hence,  a  possible  extension  of  the  methodology  presented  in  this  report  would  be  to  try  to 
caracterize  the  cost  associated  with  lack  of  reliability  from  the  resource  utilization  functions  of  each 
system  resource. 


5.4  System  design  optimization  criteria  derived  from  this 
model 

Finally,  it  has  been  described  in  Section  4  how  reliability  seems  to  be  in  contradiction  with  other 
performance  indicators.  If  it  is  assumed  that  the  performance  of  a  digital  computer  can  be 
characterized  by  a  vector,  each  component  measuring  a  different  aspect  of  performance  (for 
example,  throughput,  execution  rate,  reliability,  storage  capacity,  etc.)  the  arguments  exposed  in 
Section  4  seem  to  indicate  that  it  is  not  possible  to  raise  the  value  of  all  these  components  at  the  same 
time  (except  by  enforcement  of  fault-intolerance  in  each  resource).  Hence,  it  seems  in  principle 
possible  to  look  for  the  optimum  performance  point,  as  the  point  in  which  the  system  operates  in  a 
state  in  which  a  cost  function  associated  with  each  performance  measure  is  minimum. 
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I.  The  cyclostationary  Poisson  process 


The  problem  is  to  obtain  an  expression  for  the  pdf  of  the  time  to  failure  of  a  doubly  stochastic  failure 
process 


P(t) 


e{x«4X,TW} 


(1.1) 


where  the  time  origin  is  assumed  to  be  t  *  0  and 


\(t)  =  au(t)  +  b 

=  am(t)  +  b  +  az(t)  -  ai(t) 


(1.2) 


m(t)  being  a  periodic  function  of  time,  z(t)  a  stations  y  Gaussian  process  independent  of  m(t),  and 


2(t)  «  { 


*t> 

0  otherwise 


(13) 


Define 

A(t)  -  f  \(t)  dr 
Jo 

It  is  shown  in  [Saleh  74]  that  since 


(1.4) 


X(t)  -  ^A(t) 


(15) 


(1.1)  can  be  rewritten  as 
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(1.6) 


(1.7) 


Define 


M(t) 


f  m(t)  dt 


(18) 


Z(t) 


/  z(t)dt 


(19) 


A 


Z(t)  -  /  z(t)  dr 
Jo 


(1.10) 


The  problem  reduces  to  evaluating  the  expectation 


[aM(t)  +  aZ(t)  ■  aZ(t)  ♦  b] 


0.11) 


or,  since  m(t)  and  z(t)  are  independent,  the  problem  is  equivalent  to  evaluate 


E{e^t,M}E{et*,,a*,)1} 


(1.12) 


1.1  The  deterministic  part 

Let  us  examine  now  the  first  expectation.  m(t)  is  a  periodic  function  of  time.  However,  the  time 
origin  is,  in  principle,  unknown.  The  probability  density  function  (pdf)  given  here  is  going  to  be 
compared  with  the  estimated  pdf  of  a  real  system.  In  an  observed  pdf,  each  system  failure  will  be 
associated  with  a  time  origin  corresponding  to  the  moment  at  which  the  system  was  started.  However, 
if  the  failure  rate  is  a  time-dependent  function,  it  cannot  be  assumed  that  the  system  will  be  started 
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with  equal  probability  at  any  time  during  a  one  day  period,  in  fact,  the  system  will  be  started  more 
often  in  the  periods  of  time  in  which  it  fails  more.  Since  m(t)  is  precisely  the  periodic  component  of  the 
failure  rate, 


m 0(v)eJv 


V  4-t 
m(r)dr 


du 


(1.13) 


whehe  mQ(t)  is  m(t)  after  normalizing  to  have  area  one  in  a  one  day  period, 


mn(t) 


fflQL 


(1.14) 


m(T)dT 


and  T  is  the  period  of  m(t). 

To  evaluate  the  expectation  in  (1.13),  m(t)  will  be  approximated  by  a  finite  Fourier  series  expansion. 


m(t) 


ih  +  c0sin(nwt  +  «j>n) 


(1.15) 


where  the  following  constants  have  been  used 


(O 


1L 


cn  *  (an  +  $ 

9n  »  arctan 


(1.16) 

T 

m(t)dt 

(1.17) 

<> 1/2 
$ 

(118) 

i-ia. 

(1.19) 

bn 

T 

m(t)  cos(nwt )  dt 

(1.20) 

>„  -  /  Tn 

"  T  L 


b.  ■  /  m(t)  sln(n«t )  dt 

Jo 


(1.21) 
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Taking  into  account  this  approximation,  we  can  obtain  another  representation  for  M(t). 
that 


/ 


T+t 

m(u)di;  =  am  t  + 


EN  ac  p 

n  =  1  co$(na>T  +  <pn)  -  cos(nw(T  +  t)  +  <pn) 


] 


=  affl  t  +  ga(t,r) 


where  ga(t,r)  is  defined  as 


9a(t.T)  =  =  cosfnwr  +  9n)  -  cos(n<j(r  + 1)  +  <pn)  ] 


r  aM(t)  +  bt  -» 
e  j  can  now  be  written  as 


r  -aM(t)  +  bt  -»  (am  +  b)t  f  1  -oJLt) 

E\e  j  3  e  /  mQ(r)e  a  dT 


^a(t)e 


(am  +  b)t 


where 


/'■ 


•ga(tT) 


'a(t)  *  /  m0(r)e  a  dr 
JQ 


is  also  a  periodic  function  of  t. 


1.2  The  stochastic  part 


The  problem  is  now  to  compute 


g  r  ^a(Z(t)-Z(t)]  j 


Rrst,  note 

(1.22) 

(1.23) 

(1.24) 

(1.25) 

(1.26) 

(1.27) 


(1.28) 
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Since  z(t)  is  a  zero  mean  Gaussian  process,  aZ(t)  will  be  another  zero  mean  Gaussian  process  with 
variance  (see  [Papoulis  65],  pp.  323-325) 


<r2(t)  =  2a2 


t 


(1.29) 


where  Rzz(t1,t2)  is  the  autocorrelation  function  of  the  process  z(t).  If  in  addition  z(t)  is  stationary, 


a2(t)  =  2a2 


t 

(t-r)  R^rjdr 


(1.30) 


The  main  problem  in  the  evaluation  of  (1.28)  is  the  evaluation  of  the  statistics  of  the  process  z(t)  after 
integration,  that  is,  the  statistics  of  the  excess  area  of  a  Gaussian  process  above  a  given  level  c  in  [0,t] 
(see  Figures  1-1  and  1-2).  The  problem  of  level  crossing  for  Gaussian  processes  has  been  extensively 
treated  in  the  literature.  In  particular,  in  [Stratonovich  67]  this  problem  is  studied  in  detail,  and 
expressions  are  given  for  the  duration  of  peaks  above  a  given  level,  and  the  excess  area  under  such 
peaks.  This  is  exactly  what  is  required.  The  following  is  a  summary  of  Chapter  1-3,  Vol  II,  of 
[Stratonovich  67].  The  derivation  will  only  be  outlined  here  with  remarks  on  the  assumptions  used 
and  the  results  obtained. 


Figure  I  - 1 :  A  possible  realization  of  z(r)-z(r) 
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Figure  1-2:  A  realization  of  z(r)._The  shadowed  area  corresponds 
to  the  integral  of  z(t)  from  0  to  t,  Z(t) . 

Z(t)  is  a  random  variable  whose  exact  statistics  may  be  impossible  to  compute.  Its  value  for  a  given 
realization  of  the  process  z(t)  is  equal  to  the  addition  of  the  excess  areas  of  all  peaks  of  z(t)  above  a 
level  c.  It  is  shown  in  [Stratonovich  67],  p.  59,  that  if  the  duration  of  the  peaks  is  much  smaller  than  the 
time  between  peaks,  the  time  between  upcrossings  (downcrossings)  can  be  approximated  by  an 
exponentially  distributed  random  variable.  The  probability  of  having  k  peaks  in  [0,t]  is  then  given  by 

P(n  =  k)  =  -lajjl  e'1*  (1.31) 

where  tj  is  the  mean  number  of  peaks  per  unit  time.  In  the  case  of  a  stationary  Gaussian  process,  t|  is 
given  by  ( [Stratonovich  67],  p.  7) 

(Rj)1/2  ■£ 

e2ff2  (1.32) 


where 
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These  expressions  are  valid  only  for  "smooth"  processes.  Qualitatively,  the  smoothness  condition 
means  the  z(t)  must  be  differentiable  and  close  to  a  straight  tine  segment  during  a  sufficiently  small 
time  interval.  In  the  case  of  Gaussian  processes,  this  condition  means  that  Ru(t)  must  have  a  finite 
second  derivative  at  the  origin.  Assuming  that  this  condition  is  satisfied,  it  is  required  now  to 
characterize  the  excess  area  under  each  peak. 

If  the  duration  of  peaks  is  small  and  the  process  is  smooth,  the  second  derivative  of  z(t)  can  be 
considered  constant  over  the  duration  of  a  peak,  the  peak  can  be  assumed  to  be  of  parabolic  shape, 
and  the  excess  area  depends  only  on  the  value  of  the  first  derivative  of  z(t)  at  the  time  of  crossing  the 
level  c.  Under  these  assumptions,  the  following  expression  can  be  obtained  for  the  probability  density 
function  of  the  area  under  a  peak  ( [Stratonovich  67],  Vol  II,  pp.  68-72) 


P(s) 


f[f4 


«  vi/2,, 2/3 


(1.34) 


The  value  of  Z(t)  is  given  by 


Z(t)  *  s1  + . +  s* 


(1.35) 


where  each  of  the  s,  has  density  given  in  (1.34)  and  k  is  another  random  variable  with  density  given  in 
(1.31).  An  exact  evaluation  of  (1.28)  would  require  knowledge  of  the  joint  probability  density  function  of 
Z(t)  and  Z(t).  Since  this  impossible  to  obtain,  the  approximation  will  be  made  that 


Z(t)  =  k  E[s] 


(1.36) 


where  k  is  the  number  of  peaks  in  [0,t]  and  E[s]  is  the  mean  peak  area.  In  this  case,  Z(t)  and  Z(t)  are 
independent  random  variables  and 


r  a{Z(t)-Z(t)]  r  aZfl)',  r  -aZ(t) 

Ete  )-Eie  jE^e 


) 


(1.37) 


Since  Z(t)  is  a  zero  mean  Gaussian  variable  with  variance  <r2(t)  given  in  (1.31), 
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Et."2*} 


2v1/2a(t)  •'.00 


1 


*  am  A 

0  0  2a2(t)  dZ(t) 


(1.38) 


(139) 


and 


aim 

e 


}-E 


00  akE[s] 

kB0P<n-k>e 


(1.40) 


^(eaE(s].i), 

*  6 

p(a.c)t 

*  e 


(141) 

(1.42) 


where  p(a,c)  is  the  correction  factor 
aETal 

p(a,c)  *  i|(e  *1)  (1.43) 

Note  that  both  ij  and  E[s]  depend  on  the  value  of  the  level  c.  The  value  of  E[s]  can  be  computed  from 
(1.34). 


er  i  ,  g^rts/a 

iCfrri 

m  (2*l1/2 

4r2>,/2 

O'* 


(1.44) 


(1.45) 


If  aE[s]«<1  (that  is,  if  the  value  of  c  is  much  larger  than  the  variance  a2),  the  approximation  can  be 
made  that 
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aEfsl 

e  1 1  at  1  +  aE[s]  (1.46) 

Substituting  now  (1.46)  in  (1.46) 

p(a,c)  —  arjE[s]  (1.47) 

Substituting  now  (1.32)  and  (1.45)  into  (1.47)  a  simpler  expression  is  obtained  for  p(a,c) 

2 

2  1/2  . 

p(a,c)  *  ag  — e2a2  (1.48) 

c2 

where  a2  is  the  variance  of  z(t),  and  c  is  the  level  below  which  z(t)-z(t)  must  vanish. 


1.3  Example 

If  z(t)  is  a  zero  mean  Gaussian  process  with  variance  a2  and  autocorrelation  function 


Ra(t)  -  o  e 


0W 


(148) 


the  expectation  in  (1.40)  cannot  be  evaluated  beacuse  the  second  derivative  of  (1.49)  at  the  origin  does 
not  exist.  However,  (1.49)  can  be  approximated  by  the  following  equation 


R'a(t) 


C-Bt  -T/r.n 

e  -r^e  c] 


T/T 


2o2 


1-T^ 


(150) 


provided  that  r^<«l.  In  this  case,  the  variance  if  the  integrated  process  becomes,  according  to 
(1.30) 
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and  the  expectation  (l.i  1)  becomes 

,_2  2  -0t 

-[am  +  b  -  p(a,c)  •  [i  -  •  J 

<Aa(t)  e  &  p2 


0.51) 


(152) 


and  p(a,c)  is  given  in  equation  (1.48).  The  fact  that  the  term  (R2)1/2  cancels  in  the  approximated  value 
of  the  correction  factor  p(a,c)  has  an  interesting  physical  meaning.  R’^ft)  is  the  autocorrelation 
function  of  the  Gaussian  process  that  would  be  obtained  at  the  output  of  a  low-pass  filter  with 
bandiwth  l/rc  when  its  input  is  the  process  z(t).  The  fact  that  the  value  of  p(a,c)  is  independent  of  R2 
means  that  the  area  generated  under  the  peaks  of  z(t)  per  unit  time  is  independent  of  the  bandwith  of 
this  filter.  The  area  generated  per  peak  diminishes  for  higher  bandwith,  while  the  number  of  peaks  per 
unit  time  increases  with  the  bandwith.  Fortunately,  these  to  effects  cancel  each  other,  such  that  the 
area  generated  under  the  peaks  per  unit  time  remains  a  constant,  independent  of  the  process 
bandwith. 
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II.  Parameter  Estimation 


As  it  has  been  described  in  Section  2.5.,  the  problem  that  has  to  be  repeatedly  solved  in  this  report 
is  that  of  finding  the  minimum  of  a  function  of  the  form 


l(x )  -  23",!  H(tf(.ts,,x) - 23,1 ,  lnfh(tf.,x)] 


(H.1) 


subject  to  the  constraints 


h(tf,,x  )X)  i » 1 . ,n 


(».2) 


where  n  failures  of  a  resource  have  been  observed  at  times  tf,,  after  observing  the  system  since  ts,. 
Take,  for  instance,  the  case  of  the  distribution  obtained  from  the  simplified  cyclostationary  model 
presented  in  Section  4.1.. 


A 


...  r  r«  Ai  ) 

p(t)  «  L  o  -  y[i  -  e  J  J  e  p 


(H.3) 


Given  a  history  of  n  failures  represented  by  a  set  of  pairs  [tSj.tf.]  i « 1 . n  the  maximum  likelihood 

values  of  a,y,P  are  these  values  that  maximize  the  function 

■/{(tr.-taj 

. 1  M 


or,  equivalently,  they  are  the  values  wich  minimize  the  function 


\(a,y,P)  -  23", ,  (a-y)(tf.-ts,)  +  23". ,  -^[l  -  e  ^  ,  In  [a  -  y{1  -  e  ^ 

“  (11.5) 


subject  to  the  constraints 
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1)  «>0 

2)  fi>0 

3)  y>0 

4)  a-y>0  (11.6) 

Although  this  is  a  typical  nonlinear  programming  problem  for  which  general  methods  are  available 
(see,  for  example,  [Bazaraa  79]),  the  problem  of  minimizaing  functions  like  (11.5)  present  two 
peculiarities.  First,  every  time  that  the  function  has  to  be  evaluated  it  requires  to  compute  the  sum  of  n 
terms,  n  being  the  number  of  observed  failures.  In  the  case  of  the  file  system  described  in  section  3.1 ., 
the  number  of  observed  failures  in  five  days  of  system  operation  is  877.  Hence,  function  evaluation 
(or  gradient  evaluation)  is  computationally  expensive  and  an  efficient  method  will  be  needed. 

Second,  although  several  efficient  methods  are  known  for  minimization  subject  to  nonlinear 
inequality  constriants,  these  methods  usually  assume  that  the  constraints  are  external  to  the 
mathematical  statement  of  the  problem,  and  that  the  objective  function  can  in  fact  be  evaluated 
outside  the  constraints.  This  is  not  the  present  case.  The  fourth  constraint  in  (11.6)  must  be  satisfied 
plainly  because  the  objective  function  (11.5)  does  not  exist  unless  its  parameters  satisfy  this  constraint 
(in  the  sense  that  the  logarithm  of  a  negative  number  does  not  exist).  Indeed,  the  fourth  constraint 
says  that  the  hazard  function  must  be  positive,  and  a  solution  that  does  not  satisfy  the  fourth 
constraint  in  (11.6)  invalidates  the  existence  of  the  objective  function  itself.  Hence,  minimization 
algorithms  that  require  the  evaluation  of  the  objective  function  outside  the  constraints  cannot  be 
used. 


For  this  reason,  the  first  algorithm  to  be  used  to  find  a  minimum  of  functions  of  the  type  (11.5)  was 
the  gradient  projection  method  of  Rosen  [Rosen  60].  This  algorithm  follows  a  steepest  descent 
direction  until  one  or  several  constraints  are  violated,  projecting  then  the  gradient  on  the  subspace 
defined  by  the  active  constraints.  This  method  has  proven  to  be  very  slow  with  the  functions  tested  in 
this  report.  After  some  experimentation,  the  fastest  algorithm  found  has  been  a  slightly  modified 
version  of  a  variable  metric  algorithm  proposed  by  Powell  [Powell  78].  The  original  Powell  algorithm 
occasionally  requires  the  evaluation  of  the  objective  function  outside  the  constraints  and  has  been 
modified  such  that  the  maximum  step  size  at  each  iteration  never  leads  to  a  point  outside  the 
constraints.  The  modified  algorithm  converges  more  slowly  that  the  original  Powell  algorithm,  but  for 
ail  the  cases  in  this  report,  the  minimum  has  been  found  in  less  than  30  iterations,  which  is  a  very 
good  rate  of  convergence  given  the  functions  under  consideration. 

The  algorithm  has  been  implemented  as  an  APL  package  that  requires  the  definition  of  the 
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objective  function,  the  gradient  of  the  objective  function,  the  constraints  and  the  gradients  of  the 
constraints.  The  following  sections  desfcribe  each  function  in  detail,  providing  a  notational  dictionary 
consistent  with  the  programs  used. 
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11.1  File  System  Failures  (Cyclostationary  Model) 

Function  to  be  minimized 

KSjjk’Cdk)  *  E,.iSdk[M(ttW„  +  a 1  Ccdk  *  F^dk'^min^  "  s<flJ  ^  (tf  j-tS^) 


p  i 

P2 

■  2-".  i ln  [>*m(«r)  +  cdk ' 


-»Nw\  ] 

P 1  p2 


where 


pdk^8dk,bmln^ 


+  (w/2)1/2 


e  2(a^  ♦  otj) 


Definitions 


K,  • 

K. 

\  •  “r“i 

K»-S‘..|s 

v 

K<  ■  £|.1  K*, 

-  rrXtf,) 

(a.fi 

^•-L~ 


<a1*a-)3/2<ir/a>1/a 


e  2(0^  aj) 
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X1  *  *dk 


X2  *  Cdk 


Objective  Function 

l(x)  -  x,K,  +  [x2-x1K6-x2[-|l+-|2-]]  1C,  +  x*[  ^  ] 

* t  In  [  x,^  +  x2  -  XjKg  -  xffK^  +  Kj] 

Gradient 

dtA  *  Ki'KeK2*2xiC^+  Ka  +  + 

K5iK62xlP<31-K411 
“  “l«1  ,  O 

X1^5  +  x2  *  X1K0  *  ^3j  *  ^4|'x1 

Jl.k-E0  _ 1 _ 

°*2  xtKs  +  x2 -  XjKg  -  [Kg  +K4 ]x. 

Constraints 

0,(5?)  «Xl>0 
C2(3?)  «  x2>0 

C3(x)  «  Xlmin{K5)  +  x2  -  x,Ke  -  xfjmax^  +  K4)]  >  0 


Gradient  of  the  constraints 
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min{K,;} * 2x1  max{K_  +K.}  0 
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11.2  System  Failures  (Cyclostationary  Model) 


Function  to  be  minimized 


♦  E,”.!  [<=>», 


•  *'  sy  o2  *-  J 

P 1 


1  =  1  ^  a2  L  J 


where 


In[s3ym(tlt)  ♦  cw-p(sw,kmio)-M3swJ'0) 


*  ^SW  *  ^hw 


sy  sw 
Definitions 


-  M^-Mfts,)  K,  -  X]"m1K1( 

k2  -  2,.,^ 

Kg .  szr.,K3. 


K2f,  *  tfl*tSl 


K^.^l-e^] 


vSi-A]  k4.L>4| 


K5  -  m(tfj) 


xi  *  8» 
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X23Shw 


X3  *  Chw 


X4ak0 


Kg  as  defined  in  Section  11.1 . 


(o1+a2)3/2  (ir/2)uz 
*7  =  , 


e 


Objective  Function 


l(x)  =  (x,  +  XgjK,  +  [x3-(x1  +  x2)K6-x1K7-(x1+x2)2[-Jl+-^-]]k2 

p  1  P  2 


[  (»,**a>K5|  *  *3 - <x,  +xJ)K6-x,K7-()c,'Ht2)2tK3i  ♦  K^JJ 


Gradient 


|L  =  K,  -  KeKa  -  IV<j  -  21*,  ♦  Xj)  [  *  -g-  ]  K2  *  2|x,  ♦  Xj)  [  ^  ] 


V" 


V  VV^Py^ 


(x^+x^JK^  +  x^-ix^x^g-x^K^-lK^  +  K^Kx^+Xj) 

^  *  K,  -  K,Ka  -  2(x,  .  Xj)[  ^  ]  K,  *  2(x,  *  *,)[  i  *  Jk  ] 


V"  . 


Ks,  K6  -  -  K4,I 


(x^  +  +  x3  ‘(x^  +  X2>Kq  *  x^Ky  *  [K^  +  ](x^  ♦  Xj) 
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Jl  =  k2-EL _ 3 _ 

9x3  (x1+  x2)K5j  *  x3  -(x,  +  XgJKg  -x,^-  [K3j  +  K^Kx,  +  x^2 


|L  *  -X^fA  +  -^-]  +  £"a1 

3x4  '  T  x4  ®i  +  a2  1 


Constraints 


_ x^U/x^  +  ^Aa^qg))! 

(x1  +x2)K^  +  x3-(x1 +x2)Kg-x1K7- [Kg +K4Kx1+Xj)2 


C^x)  *  Xl>0 
C2(x)  *x2>0 
C3(x)  *  x3>0 

c4(x)  *  x4>0 

C5(x )  =  (x,  +  x2)min{K5_}  +  x2  -  (x,  +  x2)K6  -  x^  -  (x,  +  x2)2[max{K3.  +  K^}]  >  0 


Gradient  of  the  constraints 


ax,'  “  s1 . 4 


0W 

=  min{K5  }  -  V  ^  -  2xt  max{K3)  +  K^> 


-9^-  *  min^}  -  K6  -  2x ,  max^  +  K4  } 


*xiMf 


x4  q1  +  «2 


Definitions 


X,  -  «  X2  a  *1  X3  “  *1  X4  -  °2  X5  a  02 
K2  and  K2  as  in  section  11.1 . 


Objective  Function 
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£  -  ^E". ,  K*  e  ‘  A  -  ^  E?. ,  [1  -  .  ’A] 

ax  3  x3  1  ** 


*xoKj 

x  a  J  n 


+  rn  -  Z3_ _ _ _ _ 

i!l1  -x_K„  x-K, 

x^xjll-e  3  2iI-x4I1-e  5  *11 


Constraints 

C^x)  =  x,  >  0 

C2(x)  *  x,  -x2>0 
C3(x)  =  x3>0 

Gradient  of  the  constraints 


rmm 
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11.4  Simplified  Cyclostationary  Model 

Function  to  be  minimized 


l(a,Y.0)  -  E", n  (a-Y)(tfrts,)  +  E", , ^  * eP(*'**] In  [a- Y[1  - 


Definitions 

xn  ■  a  x2  ■  y  x3  *  J8 
K2  and  K2  as  in  section  11.1 . 


Objective  Function 


Kx)  *  (x1-x2)K2  +  -^E".,  [1  ln[x,  -x2[l  -e*3*2!]] 


Gradient 


JiL,,  E"  i _ 

g  =  K2-Z^u1 


Xl.x2[1-e  3  2lj 


a.  k  <t°  i,  A.r  .-."‘A 
^-Vhi-u.P 

f .  iSL,  K,  .^--aET.,  n  *  E;., 

0*3  x3  “  x. 


"*3^2 

x2K2,e  1 

x1.x2ti  e  3 


Constraints 

C,(x)  -  x,  >0 
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Gradient  of  the  constraints 


3c, 

3c, 

3x2 

3c, 

T^30 

9x3 

3c, 

3c2 

■57 

3c,  3c, 

— 2-s  0 — “-a  0 

3*3  3*i 

3c, 

3ca 

3*3 

■W.'1  J-«JI 
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11.5  Periodic  Failure  Rate 

Function  to  be  minimized 

l(sp.Cp)  -  £",1splM(tf.)-M(tsj)]  +  H". ,  Cp^-Uj)  -  H". ,  ln[spm{tf^  +  cp] 


Definitions 


xi  “8p 


X2*Cp 


Kr  K2,  and  Kg  as  defined  in  Section  11.1 


Objective  Function 


l(x)  -  x,K,  +  x2K2-X]"s1  ln[XlKj.  +  xj 


Gradient 


K- 

JL  »  K,  -  S,  . - ^L 

3x1  X1K5,+ 


5|  +  *2 


»  K2-£?  . - 1— 

2  x.K-  +  x. 


0*2 


5j  T  2 


Constraints 


parameter  estimation 


Gradient  of  the  Constraints 
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III.  Autocorrelation  function  estimation 


Given  an  ergodic  and  stationary  process  z(t),  the  problem  is  to  estimate  the  function 


lim 

Tt00  2T 


z(t  +  r)z(t)  dt 


(111.1) 


For  a  finite  record  of  observed  values  z(n),  the  autocorrelation  function  is  usually  estimated  using  the 
expression 


Fyn) 


1  VN-n 

-J-2wi  =  1z(i  +  n)z(i) 


(IH*2) 


This  estimate  is  intuitive  except  for  the  factor  1/n.  Since  N-n  terms  are  summed,  it  seems  that  1  /(N-n) 
would  be  more  exact.  In  fact  (111.2)  is  a  biased  estimator  of  the  real  autocorrelation  function.  However, 
its  expected  error  is  smaller  than  the  expected  error  that  would  be  obtained  using  the  (unbiased) 
estimator  with  factor  l/(N-n)  [Jenkins  68]. 

In  the  cases  presented  in  this  report  the  values  of  z(n)  are  not  directly  observable.  In  the  case  of 
sampling  the  values  of  fraction  of  time  in  Kernel  mode,  what  was  measured  was  the  average  fraction 
of  time  in  Kernel  mode  during  the  last  second,  recording  a  sample  every  five  minutes.  In  the  case  of 
the  number  of  blocks  accessed  to  the  file  system,  the  available  samples  are  the  number  of  blocks 
accessed  during  the  last  five  minutes,  also  with  a  resolution  of  five  minutes.  The  measured  values  are 
not  the  values  of  z(n),  but  the  values  of  the  process 


z‘(n) 


z(t)  dt 


where  A  equals  five  minutes  or  one  second,  and  the  available  samples  of  z’(n)  are  five  minutes 
appart. 


It  has  been  observed  that  in  the  two  cases  studied  in  this  report,  the  autocorrelation  function 
suggests  an  approximation  of  the  form 


■pM  pJH 
Ra(t)  ■  a,e  1  +  o2e  * 


(111*4) 


The  problem  is  then  to  estimate  the  values  of  the  av  P{  from  the  observed  values  of  z'(n).  If 
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Rz,z,(t)  -  «’/lW  +  a'2e^ 


it  is  easy  to  show  that 


RH(t)  *  o^e  +  o2e  ? 


(Hl*5) 


(111*6) 


where 


_ 4!S_ 

'  2(cosh(0.  A)  - 1] 


(IH*7) 


The  problem  is  then  to  estimate  the  values  of  the  a’.,  0.  using  (111.2)  and  the  observed  values  of  z’(n), 
and  use  (111.7)  to  obtain  the  values  of  a,  of  the  autocorrelation  function  of  z(t). 

Unfortunately  it  has  not  been  possible  to  follow  this  procedure.  The  accuracy  of  the  estimated 
autocorrelation  function  is  limited  basicly  by  two  factors  :  the  sampling  frequency  and  the  length  of 
the  available  record,  N.  Although  many  techniques  exist  for  power  spectrum  estimation  that  take  into 
account  these  two  factors  [Oppenheim  75]  (the  power  spectrum  is  the  Fourier  transform  of  the 
autocorrelation  function),  no  techniques  are  available  for  correcting  the  estimates  of  the 
autocorrelation  function  itself. 

If  the  sampling  frequency  is  comparable  to  the  bandwith  of  the  power  spectrum,  the  power 
spectrum  estimate  may  be  poor  due  to  aliasing.  Under  these  conditions,  the  estimate  of  the 
autocorrelation  function  given  by  (111.2)  may  take  negative  values.  This  is  precissely  what  happens  for 
the  estimated  autocorre'.ation  function  of  the  file  system  utilization  process  as  shown  in  Figure  3-2. 
For  a  sampling  frequency  equal  to  one,  the  bandwith  of  the  process  would  be  equal  to  =  0.59,  that 
is,  the  sampling  frequency  is  not  even  twice  the  process  bandwith. 

The  solution  adopted  has  been  to  estimate  the  a.  and  0.  directly  from  a  history  of  failures  as 
described  in  Section  4.1 .,  and  to  estimate  the  variance  of  the  process  a2  .  Since 


a 


2 

zz 


*  +  a 


2 


(IH*8) 


and 
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And  knowing  ,  p.  and  <t^(  the  values  of  the  can  be  computed 
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